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The  proposed  theory  of  causally  structured  discrete  fields  studies  integer  values  on  directed  edges 
of  a  self-similar  graph  with  a  propagation  rule,  which  we  define  as  a  set  of  valid  combinations  of  inte¬ 
ger  values  and  edge  directions  around  any  vertex  of  the  graph.  There  is  an  infinite  countable  number 
of  variants  of  the  theory  for  a  given  self-similar  graph  depending  on  the  choice  of  propagation  rules, 
some  of  these  models  can  generate  infinite  uncountable  sets  of  patterns.  This  theory  takes  minimum 
assumptions  of  causality,  discreteness,  locality,  and  determinism,  as  well  as  fundamental  symmetries 
of  isotropy,  CPT  invariance,  and  charge  conservation.  It  combines  the  elements  of  cellular  automata, 
causal  sets,  loop  quantum  gravity,  and  causal  dynamical  triangulations  to  become  an  excellent  can¬ 
didate  to  describe  quantum  gravity  at  the  Planck  scale.  In  addition  to  the  self-consistent  generation 
of  spacetime  and  metrics  to  describe  gravity  and  an  expanding  closed  Universe,  the  theory  allows  for 
the  many-worlds  interpretation  of  quantum  mechanics.  We  also  demonstrate  how  to  get  to  unitary 
evolution  in  Hilbert  space  for  a  stationary  Universe  with  deterministic  propagation. 

PACS  numbers:  04.60.-m,  04.60.Nc,  04.90.+e,  03.65.-w,  03.65.Ta,  03.65. Ud 


I.  INTRODUCTION 

There  is  a  long-standing  idea  that  the  Universe  must  be 
deterministic,  discrete,  and  local  at  the  Planck  scale  T], 
There  is  also  a  temptation  to  build  such  a  Universe  via 
cellular  automata  that  can  demonstrate  complex  behav¬ 
ior  even  from  simple  rules  @| ,  but  the  patterns  from  cellu¬ 
lar  automata  have  a  global  time  axis,  which  is  not  appli¬ 
cable  to  describe  general  relativity.  Another  idea  is  that 
the  Universe  may  be  a  tiling  of  spacetime  Q.  We  intro¬ 
duce  a  theory  that  combines  these  two  ideas  into  causal 
tilings  of  spacetime  that  represent  patterns  of  integer  val¬ 
ues  on  edges  of  self-similar  directed  acyclic  graphs. 

On  a  historical  note,  already  at  the  time  of  the  foun¬ 
dation  of  quantum  mechanics  and  general  relativity, 
there  was  a  suggestion  that  the  Universe  is  a  compu¬ 
tation  0.  Later,  several  authors  discussed  this  idea  fur¬ 
ther  Q ,  @ ,  [?] .  The  causal  set  theory  was  the  first  to  use 
the  idea  of  causality  to  study  the  structure  of  spacetime 
at  the  Planck  scale  [§j,  but  it  required  external  dynamics 
to  reproduce  physically  meaningful  results.  Casual  dy¬ 
namical  triangulation  allows  to  create  physically  mean¬ 
ingful  causal  tilings  of  spacetime  Q,  but  it  required  an 
assumption  for  the  Hilbert-Einstein  action  with  model 
parameters  a  and  A.  One  of  the  closest  predecessors  to 
the  proposed  theory  is  loop  quantum  gravity  and  its  rep¬ 
resentation  with  spin  newtworks  0  .  |ll|.  However,  this 
theory  again  required  external  assumptions  and  started 
from  a  formulation  of  general  relativity  with  new  param¬ 
eters.  The  approach  of  this  paper  is  to  start  purely  from 
discrete  geometry  and  to  search  for  ways  how  general 
relativity  and  quantum  mechanics  emerge  from  there. 

We  organize  the  paper  in  the  following  way.  First, 
we  introduce  a  path-integral  approach  for  a  general  case 
of  discrete  integer  fields  on  edges  of  a  self-similar  graph 
(lattice)  in  any  dimension.  We  postulate  discrete  fields 
as  integer  values  and  edge  directions  of  the  graph.  We  in¬ 
troduce  a  propagation  rule  that  describes  all  valid  combi¬ 
nations  of  fields  on  edges  around  any  vertex  of  the  graph. 


We  allow  only  for  a  single  initial  vertex. 

Second,  we  introduce  a  toy  model  in  two-dimensions 
on  a  honeycomb  graph  to  show  the  intuition  behind  the 
theory.  We  use  a  Matlab  code  that  generates  and  visu¬ 
alizes  two-dimensional  patterns  on  this  graph  (12  .  We 
show  examples  of  patterns  that  are  generated  by  this  code 
throughout  the  paper  in  order  to  support  the  discussion. 

Third,  we  discuss  the  phenomenology  that  can  be 
demonstrated  by  this  theory.  We  describe  a  family  of  gen¬ 
eralized  diamond  crystal  lattices  that  allow  for  isotropic 
pattern  generation,  which  is  crucial  for  the  representa¬ 
tion  of  a  uniform  space.  We  describe  the  Big  Bang  that 
originates  causally  from  an  initial  vertex.  We  demon¬ 
strate  the  existence  of  black  holes  at  the  Planck  scale. 
We  introduce  the  future  and  past  light  cones.  We  pro¬ 
pose  a  covariant  definition  of  a  quantized  spacetime  in¬ 
terval  purely  from  geometry.  This  measure  adds  scale  to 
the  causal  relationship  between  vertices  of  the  graph  and 
therefore  allows  to  reconstruct  a  continuous  Riemannian 
manifold.  We  postulate  the  energy-momentum  tensor 
together  with  the  cosmological  term  A  from  the  Einstein 
field  equations.  We  introduce  CPT  invariance  and  charge 
conservation  within  the  theory.  We  make  a  hypothesis 
that  elementary  particles  and  low  energy  quantum  field 
theories  may  arise  from  the  dynamics  of  the  patterns  on 
the  graph.  We  suggest  that  non-computable  charge  con¬ 
serving  propagation  rules  may  be  relevant  for  the  creation 
of  mass  and  existence  of  inertia.  We  demonstrate  the 
way  to  calculate  entropy  for  a  case  of  a  stationary  Uni¬ 
verse.  We  briefly  discuss  the  relevance  of  this  theory  to 
describe  quantum  phenomena  such  as  Heisenberg  uncer¬ 
tainty,  the  double-slit  experiment,  EPR  paradox,  Bell’s 
theorem,  and  quantum  entanglement.  For  this  matter, 
we  refer  to  the  many-worlds  interpretation  of  quantum 
mechanics.  Finally,  we  show  a  way  to  get  to  unitary 
evolution  in  a  Hilbert  space  in  the  case  of  a  stationary 
Universe  with  deterministic  propagation. 

We  conclude  with  a  discussion  of  possible  ways  to  test 
the  theory  in  order  to  make  it  falsifiable. 
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II.  GENERAL  THEORY  WITH 
PATH-INTEGRAL  FORMULATION 

The  proposed  theory  describes  all  valid  configura¬ 
tions  of  integer  values  and  edge  directions  on  self-similar 
graphs  (lattices)  given  a  propagation  rule.  We  start  with 
important  definitions: 

Definition  1  A  translation  of  a  graph  is  a  permutation 
of  its  vertices  that  preserves  connectedness  between  all 
vertices  of  the  graph. 

Definition  2  A  self-similar  graph  has  all  vertices  with 
k  direct  neighbors,  and  it  allows  for  a  translation  that 
places  any  vertex  i  instead  of  any  other  vertex  j. 

Definition  3  Any  edge  of  a  self-similar  graph  has  a  class 
(angular  direction)  assigned  to  it:  c  =  1  so  that 

any  translation  does  not  mix  classes  of  edges  and  moves 
all  edges  of  one  class  to  another  class. 

An  example  of  a  self-similar  graph  is  a  cubic  lattice  in 
n-dimensions.  There  are  k  =  2d  direct  neighbors  for  each 
vertex  and  c  =  2d  classes  of  edges  two  for  each  direction 
along  each  dimension  d  =  n.  Cubic  lattices  are  simple 
and  easy  for  modeling,  but  they  may  be  irrelevant  for  not 
having  enough  isotropy. 

Another  example  is  a  general  diamond  crystal  lattice  in 
n-dimensions,  which  is  a  subset  of  vertices  of  an  (n+1)- 
dimensional  cubic  lattice  with  integer  coordinates  that 
sum  either  to  one  or  zero.  The  edges  for  this  graph 
are  formed  between  direct  neighbors  in  the  initial  (n+1)- 
dimensional  cubic  lattice.  A  general  diamond  crystal  lat¬ 
tice  in  d  =  n  dimensions  has  k  =  d  +  1  direct  neighbors 
for  each  vertex  and  c  =  d  + 1  classes  of  edges  one  for  each 
direction  of  edges.  A  general  diamond  crystal  lattice  is 
a  honeycomb  graph  in  two  dimensions  and  a  diamond 
crystal  lattice  in  three  dimensions. 

Definition  4  An  integer  field  variable  ipij  = 
— m, . . . ,  —  1, 0, 1, . . . ,  m  is  defined  for  any  two  adja¬ 
cent  vertices  i  and  j,  where  a  positive  integer  value 
means  direction  from  i  to  j,  and  a  negative  integer 
value  means  direction  from  j  to  i.  These  fields  are 
antisymmetric:  ipij  =  —ipji- 

Zero  value  of  ipij  means  that  the  edge  does  not  have 
an  assigned  integer  value  and  direction.  Some  variations 
of  the  theory  may  forbid  zero  integers  on  edges.  In  this 
case,  all  edges  will  have  an  assigned  direction. 

Definition  5  Any  vertex  i  has  a  set  of  interger  fields 
on  its  edges  that  are  defined  by  a  tuple  Ti  = 
{^i,vU,i)i  ■  ■  ■  i  V’iMi.fc)};  with  a  vicinity  function  j  = 
v(i,  c )  that  orders  edges  ipij  by  classes  c  =  1, . . . ,  k. 

Definition  6  A  propagation  rule  V  is  defined  for  any 
vertex  i  as  a  subset  of  all  valid  tuples  of  integer  fields  on 
edges:  V  C  V^i, „(*,!)>  •  •  • ,  ipi,v(i,k)}- 


Definition  7  A  vertex  i  satisfies  a  propagation  rule  V 
if  the  tuple  of  integer  fields  around  the  vertex  7)  belongs 
to  the  subset  of  all  allowed  tuples  of  integer  fields  for  the 
propagation  rule:  71  GV. 

Definition  8  We  introduce  three  causality  constraints 
for  the  discrete  fields  ipij : 

•  Propagation  rule:  Vi  :  7i  €  V 

•  No  cycles :  $ii,...,in  :  ipili2  >  0, . . . ,  ipinil  >  0 

•  Unique  initial  vertex:  3!6  :  Vj  ipbj  >  0 

These  three  constraints  make  sure  that  the  edges  of 
the  graph  form  a  directed  acyclic  graph  (DAG)  and  all 
vertices  of  the  graph  satisfy  a  specific  propagation  rule 
V. 

Now  we  can  formulate  the  main  idea  of  the  proposed 
theory: 

The  casual  discrete  field  theory  (CDFT)  stud¬ 
ies  properties  of  integer  value  fields  ipij  on 
edges  of  self-similar  graphs  under  causality 
constraints  with  a  given  propagation  rule  V . 

We  define  any  measurable  physical  quantity  as  a  func¬ 
tional  over  field  configurations:  F[ip],  Also  for  conve¬ 
nience,  we  define  a  notation  for  a  constraint  functional 
Cv,b\jP\  for  a  certain  propagation  rule  V  and  an  initial 
vertex  b.  This  constraint  functional  is  proportional  to 
the  likelihood  of  each  configuration  of  fields  that  satisfy 
causality  constraints  and  0  otherwise.  This  constraint 
functional  C-p,b  [ip\  is  a  discrete  analog  of  the  action  term 
gtSM  and  time  ordering  for  path  integrals  in  quantum 
field  theories. 

We  define  an  expectation  value  of  an  operator  F  with 
an  analog  of  path  integral  summation  over  all  valid  field 
configurations  undger  causality  constraints: 

{F)  ~  TWcUW]  (1) 

The  theory  is  constructed  with  self-similar  graphs  be¬ 
cause  they  are  composed  of  vertices  with  identical  local 
topology.  Any  vertex  can  be  replaced  by  any  other  vertex 
via  translations  so  that  the  graph  would  entirely  repeat 
itself.  The  theory  can  be  easily  generalized  even  further 
by  allowing  lattices  with  several  classes  of  vertices  with 
different  local  topology  and  different  propagation  rules 
for  each  class,  but  we  keep  it  out  of  the  discussion  for 
simplicity. 

A  propagation  rule  for  a  given  self-similar  graph  is  the 
cornerstone  of  the  proposed  theory  because  it  fully  de¬ 
fines  all  patterns  of  integer  fields  and  edge  directions.  A 
propagation  rule  is  a  discrete  analog  of  interactions  be¬ 
tween  fields  in  continuous  field  theories. 

Definition  9  Any  propagation  rule  V  for  a  self-similar 
graph  with  k  nearest  neighbors  for  each  vertex  is  encoded 
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by  a  pair  of  numbers  m  :  E,  where  m  is  a  maximum  inte¬ 
ger  value  on  edges  and  E  is  an  encoding  of  a  subset  of  all 
numbers  in  (2m+l)-base  numeric  system  with  k  digits. 
Thus  the  encoding  E  is  a  number  from  0  to  2(2m+1)  —  1. 

Encoding  of  all  propagation  rules  with  two  positive 
integer  numbers  means  that  there  is  an  infinite  countable 
number  of  different  models  for  a  given  self-similar  graph. 

We  use  integer  values  on  edges  to  keep  the  theory  dis¬ 
crete  with  countable  variables,  but  the  most  general  at 
the  same  time.  We  use  a  limit  of  m  for  these  values  to 
avoid  infinite  parameter  space  for  propagation  rules.  A 
more  general  theory  is  possible  with  infinite  countable  or 
continuous  sets  of  values  on  edges  with  some  functional 
structure  for  propagation  rules,  but  this  goes  beyond  the 
scope  of  this  paper. 

We  have  defined  the  most  general  propagation  rule 
V .  This  propagation  rule  may  have  additional  limita¬ 
tions  and  symmetries  like  the  absence  of  zero  value  edges, 
isotropy,  CPT  invariance,  and  charge  conservation.  We 
consider  these  constraints  in  more  detail  in  the  discussion 
section. 

A  specific  propagation  rule  V  for  a  given  self-similar 
graph  defines  a  CDFT  model.  Practical  calculations 
show  that  even  models  with  small  m  can  generate  com¬ 
plex  causally  related  patterns  or,  in  other  words,  tilings 
of  space. 

Definition  10  Each  valid  configuration  of  integer  fields 
on  edges  corresponds  to  a  tiling  of  space,  where  each  tile 
corresponds  to  a  vertex  with  a  particular  tuple  71  from 
valid  tuples  in  the  propagation  rule  V .  In  this  tiling  two 
tiles  Ti  and  Tj  may  be  adjacent  only  if  they  share  an  edge 
with  the  same  integer  value  ipij  =  —ipji- 

From  the  perspective  of  tilings,  a  propagation  rule  is 
a  set  of  valid  tiles  that  can  tile  the  whole  space,  and  any 
two  adjacent  tiles  have  to  have  matching  absolute  values 
of  positive  and  negative  integers  from  1  to  m. 

The  presence  of  the  initial  vertex  b  is  essential  to  define 
a  starting  point  for  patterns  on  the  graph,  whereas  mul¬ 
tiple  initial  vertices  will  add  too  many  degrees  of  freedom 
in  this  theory.  There  is  also  a  sequential  growth  interpre¬ 
tation  for  patterns  on  the  graph.  Each  field  configuration 
on  a  graph  allows  at  least  one  sequential  enumeration  of 
vertices  so  that  each  new  vertex  has  all  its  directed  edges 
from  already  enumerated  adjacent  vertices. 

Some  propagation  rules  allow  for  the  computation  of 
patterns  by  sequential  growth  of  tiles  one  by  one  with  all 
valid  local  combinations.  These  propagation  rules  may 
generate  only  periodic  patterns  or  may  allow  for  more 
complex  patterns.  Also,  there  are  remaining  propaga¬ 
tion  rules  that  forbid  certain  combinations  of  values  on 
incoming  edges.  In  this  situation,  a  sequential  placement 
of  tiles  may  end  up  with  a  stalemate  when  there  is  no 
valid  tile  for  the  next  vertex.  The  propagation  rules  of 
the  second  type  may  not  have  a  possible  tiling  at  all, 
or  they  may  generate  non-computable  tilings  in  a  sense 
defined  by  Roger  Penrose  Q,  who  had  ideas  related  to 


the  proposed  theory,  but  his  famous  tilings  lacked  causal 
structure  which  is  essential  for  the  CDFT. 

Definition  11  A  propagation  rule  is  trivial  if  it  gener¬ 
ates  only  periodic  tilings  or  does  not  generate  any  tiling 
at  all. 

Definition  12  A  propagation  rule  is  computable  if  it  al¬ 
lows  for  a  deterministic  construction  of  tilings  by  the  se¬ 
quential  placing  of  tiles  one  by  one  in  a  recursion  by  ver¬ 
tices  starting  from  the  initial  vertex  b  and  by  choosing 
from  all  valid  tiles  for  each  next  vertex. 

Non-trivial  computational  propagation  rules  are  useful 
for  studying  because  they  allow  for  the  computational 
generation  of  patterns  and  calculation  of  desired  metrics 
and  properties.  However,  it  is  highly  likely  that  non¬ 
trivial  non-computational  propagation  rules  are  more  rel¬ 
evant  for  a  description  of  physical  reality,  because  they 
intrinsically  non-local  globally:  a  local  change  in  a  type 
of  tile  for  a  vertex  i  may  put  constraints  on  a  type  of  tile 
for  another  vertex  j  that  does  not  have  a  direct  causal 
relationship  with  the  vertex  i.  This  situation  resembles 
quantum  non-locality,  which  we  cover  in  more  detail  in 
the  discussion  section. 

All  possible  patterns  for  a  given  graph  and  a  given 
propagation  rule  form  a  set  of  overlapping  patterns.  Each 
pair  of  patterns  in  this  set  have  “joint  past”,  which  is 
a  maximum  number  of  connected  vertices  identical  for 
both  patterns.  As  a  result,  this  set  defines  all  possible 
evolutions  for  a  given  variant  of  a  causal  discrete  field 
theory.  We  call  this  set  of  patterns  a  Multiverse,  which 
contains  all  copies  of  possible  Universes.  We  postulate 
that  each  copy  of  a  Universe  participates  in  path  integral 
summation  over  field  configurations  in  Eq.  [l] 

There  are  propagation  rules  that  can  generate  an  in¬ 
finite  uncountable  number  of  Universes  within  a  Multi- 
verse  covering  a  continuous  spectrum  of  possibilities.  We 
can  demonstrate  that  for  a  simple  propagation  rule  with 
to  =  2  and  all  possible  incoming  and  outgoing  directions 
of  edges.  Even  for  each  fixed  directed  acyclic  graph,  there 
is  an  uncountable  number  of  patterns  of  integer  values  of 
1  or  2  on  edges,  because  this  corresponds  to  all  possible 
subsets  of  an  infinite  countable  set  of  directed  edges. 

The  proposed  CDFT  program  for  future  research  is  to 
study  the  properties  of  non-trivial  patterns  for  different 
self-similar  graphs  and  propagation  rules.  In  the  remain¬ 
ing  part  of  this  paper,  we  discuss  why  CDFT  models  are 
relevant  for  the  description  of  physical  reality. 


III.  A  TWO-DIMENSIONAL  EXAMPLE  FOR 
ILLUSTRATION 

We  use  two  simple  models  in  order  to  visualize  the  re¬ 
sults  of  the  discussion.  We  take  a  hexagonal  honeycomb 
graph  in  two  dimensions.  The  first  model  is  the  most 
simple  one.  It  has  only  one  integer  1  on  the  edges  of  the 
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FIG.  1.  An  example  of  a  propagation  rule  for  two  simple 
models  on  a  honeycomb  graph  in  two  dimensions.  Model  A: 
four  valid  types  of  vertices  in  sector  A.  Model  A-D:  sixteen 
valid  types  of  vertices  in  sectors  A-D.  Each  type  of  vertex 
also  allows  all  possible  rotations.  Edges  with  the  same  integer 
value  have  the  same  color:  1  -  black,  2  -  red,  3  -  blue. 


graph  and  all  four  valid  types  of  vertices  with  all  possi¬ 
ble  directions  of  edges  including  all  rotations,  see  Model 
A  on  Figure  [T]  The  second  model  is  a  little  bit  more 
complicated.  It  has  three  integers  1  —  3  on  the  edges  of 
the  graph  and  sixteen  valid  types  of  vertices  including  all 
rotations,  see  Model  A-D  on  Figure  [H 

An  expample  of  enconding  m  :  E  for  the  simplest  pos¬ 
sible  Model  A  is  1  :  85197125,  where  the  second  number 
equals  to  2®,  where  i  £  [0,2,2*3,2*3  +  2,2*9,2*9  + 
2, 2*9  +  2*3,2*9  +  2*3  +  2]. 

Model  A  is  non-trivial  and  computable  because  it  al¬ 
lows  for  all  possible  directions  for  a  directed  acyclic 
graph.  We  use  a  Matlab  code  to  generate  pattern  evolu¬ 
tion  for  Model  A  [T^].  We  start  from  the  initial  vertex 
b  and  sequentially  add  vertices  with  valid  types  one  by 
one.  We  pick  each  next  vertex  randomly  from  all  possi¬ 
ble  locations  and  valid  vertex  types  with  equal  likelihood. 
Figure  [2]  shows  an  example  for  several  first  steps  of  evolu¬ 
tion  for  Model  A  with  vertices  not  far  than  50  links  from 
the  initial  vertex  b. 

Model  A-D  is  also  non-trivial  because  it  contains  all 
valid  vertices  of  Model  A.  Model  A-D  is  non-computable 
because  it  does  not  allow  all  possible  combinations  for 
incoming  edges.  The  reason  for  having  Model  A-D  as  an 
additional  example  will  be  apparent  during  the  discussion 
of  symmetries  for  propagation  rules. 


IV.  DISCUSSION 

Planck  element  of  spacetime.  We  now  switch  from 
pure  mathematics  to  the  discussion  of  physical  phenom¬ 
ena  that  may  be  relevant  for  the  proposed  theory.  The 
essential  foundation  of  it  is  the  concept  of  causality  and 
discreteness  of  spacetime.  We  start  with  several  essential 


FIG.  2.  A  pattern  on  a  honeycomb  graph  with  Model  A. 
The  chart  shows  the  result  of  propagation  when  the  first  tile 
have  reached  a  vertex  50  edges  away  from  the  initial  vertex  b. 
Each  triangle  tile  corresponds  to  the  nearest  neighbourhood 
of  a  vertex.  Colors  show  the  order  of  propagation. 

definitions. 

Definition  13  Each  vertex  v  of  the  graph  is  postulated 
to  be  an  atomic  Planck  element  of  spacetime.  Later  on, 
these  vertices  are  called  spacetime  events. 

Definition  14  A  chain  of  links  l(v,w )  between  two 
events  v  and  w  is  any  connected  sequence  of  edges  of 
the  graph  regardless  of  their  direction. 

Definition  15  A  path  of  links  l±(v,w)  between  two 
events  v  and  w  is  a  chain  with  all  its  edges  directed  from 
v  to  w  for  l+(v,  w)  or  vice  versa  from  w  to  v  for  l-  (v,  w). 

Definition  16  The  length  of  a  chain  |/(u,w)|  or  a  path 
|/±(u,u))|  is  defined  as  a  number  of  links  in  it. 

Definition  17  Every  vertex  v  has  its  past  w  £  J~(v )  : 
3l-(v,w)  and  its  future  w  £  J+(v)  :  3 l+(v,w). 

These  two  sets  of  events  are  analogs  of  light  cones  orig¬ 
inating  from  an  event  v  into  the  past  and  future,  see 
Figure  [3] 

The  size  of  the  Universe.  For  each  event  v  we 
define  a  current  size  of  the  Universe,  which  does  not  have 
physical  meaning,  but  essential  for  the  following  analysis. 

Definition  18  A  cosmological  epoch  is  a  set  of  events 
with  equal  minimal  paths  from  the  initial  event  b:  w  £ 
E{n )  :  min(|Z_|_(&,  ie)|)  =  n. 

By  definition,  any  two  events  from  the  same  epoch  have 
at  least  the  initial  event  &  as  a  joint  predecessor.  The 
number  of  events  in  a  certain  epoch  defines  the  volume 
of  the  Universe  at  this  stage.  See  Figure[3]for  an  example 
of  two  epochs  at  n  =  1500  and  n  =  2000  for  Model  A. 
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FIG.  3.  An  example  of  two  epoches  with  n  =  1500  and 
n  =  2000  (circle  tiles)  with  future  and  past  light  cones  start¬ 
ing  from  an  event  from  the  first  ephoch  (triangle  tiles).  An 
evolution  is  shown  for  Model  A  only  for  events  with  causal 
descendants  at  epoch  n  =  2500  (outside  of  interiors  of  black 
holes) . 

The  Big  Bang  naturally  occurs  for  non-trivial  CDFT 
models  as  events  of  the  future  epochs  expand  from  the 
initial  event  b.  There  is  no  global  time  in  the  model, 
which  is  a  desired  property  for  a  covariant  theory  of  gen¬ 
eral  relativity. 

The  same  propagation  rule  defines  the  vicinity  of  any 
event,  so  the  Big  Bang  is  happening  at  every  event, 
which  follows  the  ideas  of  Andrei  Linde  on  the  inflation¬ 
ary  model  of  the  Universe  (l3|.  The  proposed  theory 
describes  an  eternally  expanding  closed  Universe  with  a 
possibility  that  some  regions  collapse  into  giant  gravita¬ 
tional  singularities,  so  that  local  observers  may  perceive 
them  as  closed  Universes  with  a  finite  lifetime. 

Black  holes.  Some  of  the  events  may  not  have  causal 
descendants  at  all.  Suppose  a  vertex  i  has  only  incoming 
links  V)  :  ipij  <  0.  Some  of  the  causal  predecessors  of 
this  event  i  also  do  not  have  descendants  further,  so  all 
these  events  is  an  analog  of  an  interior  of  a  black  hole  at 
the  Planck  scale.  The  event  horizon  for  such  a  black  hole 
is  a  set  of  links  that  separate  events  without  descendants 
further  than  the  event  i  from  events  with  descendants 
further  than  the  event  i. 

The  event  horizon  is  not  a  special  object  in  this  theory, 
and  there  is  no  a  “firewall” .  The  paths  of  links  that  cross 
this  border  do  not  experience  any  change  in  local  dynam¬ 
ics.  The  only  consequence  is  that  these  paths  will  never 
reach  infinity  and  they  will  collapse  into  a  singularity  at 
the  center  of  a  black  hole  after  a  finite  number  of  steps. 

The  proposed  theory  generates  many  singularity 
events,  see  Figure[4]for  example.  Spacetime  at  the  Planck 
scale  forms  quantum  foam  with  high  fluctuations  of  phys¬ 
ical  quantities  and  high  density  of  the  Planck  scale  black 
holes.  In  this  example,  propagation  that  reaches  infinity 


FIG.  4.  A  demonstration  of  black  hole  singularities.  Evo¬ 
lution  is  shown  for  Model  A.  Propagation  is  shown  for  200 
epochs  for  events  (triangle  tiles)  with  causal  descendents  at 
epoch  2500  (outside  of  interiors  of  black  holes).  Circles  show 
singularity  events  with  no  causal  descendents. 

concentrates  within  a  few  narrow  bands  of  events  that 
are  separated  by  pockets  of  interiors  of  large-scale  black 
holes.  In  four  dimensions,  the  events  that  reach  infinity 
are  separated  by  four-dimensional  pockets  correspond¬ 
ingly.  The  cosmological  inflation  of  the  Universe  should 
keep  these  narrow  bands  long  enough  for  the  size  of  the 
Universe  to  grow  much  faster  than  it  expands. 

Isotropy  of  space.  The  real  Universe  demonstrates 
the  isotropy  of  space  so  that  all  directions  are  equivalent. 
Any  good  theory  on  a  lattice  has  to  demonstrate  this 
property  too.  There  are  for  sure  graphs  and  propagation 
rules  in  the  proposed  theory  that  do  not  have  isotropic 
expansion  from  the  initial  event  b.  For  example,  in  the 
case  of  a  simple  cubic  lattice  in  n-dimensions,  expansion 
with  any  propagation  rule  seems  to  preserve  the  shape  of 
a  cube  with  softened  edges. 

General  diamond  crystal  lattices  are  more  complex 
than,  e.g.,  n-dimensional  cubic  lattices,  but  they  are 
highly  isotropic  and  have  closer  symmetry  to  a  uniform 
space  0.  These  lattices  are  unique  because  they  have 
the  minimum  number  of  edges  for  each  vertex,  and  these 
edges  stretch  uniformly  in  an  n-dimensional  space  that 
embeds  the  lattice.  Also,  there  is  a  translation  of  the 
graph  that  can  exchange  every  two  edges  from  the  same 
vertex,  which  is  not  the  case  for  a  cubic  lattice.  Every 
vertex  in  a  general  diamond  crystal  lattice  is  a  center 
of  an  n-dimensional  simplex.  The  honeycomb  graph  is  a 
two-dimensional  case  of  a  general  diamond  crystal  lattice 
with  a  triangle  simplex. 

The  propagation  rule  also  has  to  be  isotropic,  so  that 
local  topology  of  edges  around  vertices  does  not  depend 
on  the  direction  of  a  class  of  edges. 

Definition  19  A  propagation  rule  V  is  isotropic,  if  for 
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each  tuple  71  =  {ipi  v(i, 1),  ■  ■  ■ ,  ^  contains  all  its 

rotations  with  respect  to  directions  for  each  edge  class. 
For  a  general  diamond  lattice  the  rotations  of  a  tuple 
mean  all  parity  preserving  permutations  of  edge  classes 
within  the  tuple. 

There  are  at  least  two  ways  to  check  isotropy  of  pat¬ 
terns  generated  by  a  CDFT  model.  The  first  approach 
is  to  check  numerically  that  the  fluctuations  of  Sr/r  are 
vanishing  for  epochs  with  large  n.  Here  we  take  r  as  a 
Euclidean  distance  between  a  vertex  of  an  epoch  and  the 
initial  vertex  b  in  an  n-dimensional  Euclidean  space  that 
embeds  the  lattice.  We  mean  that  a  lattice  is  embedded 
into  a  Euclidean  space  of  a  minimum  dimension  when  the 
distance  between  adjacent  vertices  is  equal  to  one,  and 
the  angles  between  any  two  pairs  of  edges  are  the  same  if 
there  is  a  translation  of  the  graph  that  exchanges  these 
pairs  of  edges.  The  second  approach  purely  follows  from 
the  structure  of  the  graph.  One  can  derive  spacetime 
metrics  from  a  directed  acyclic  graph  of  links  and  check 
its  isotropic  property  on  a  large  scale.  Both  approaches 
are  computationally  demanding  and  yet  to  be  explored 
for  high  dimensions. 

Lorentz  invariance.  The  link  between  Lorentzian 
spacetime  manifolds  and  causal  sets  is  known  for  quite 
some  time.  Spacetime  events  in  a  Lorentzian  spacetime 
manifold  form  a  causal  set  |l5j.  In  the  opposite  direc¬ 
tion,  the  mainstream  causal  set  theory  (CST)  has  a  cen¬ 
tral  conjecture  (also  called  the  Hauptvermutung)  that  a 
Lorentzian  manifold  can  approximate  any  causal  set  with 
specific  properties  in  a  unique  way  under  isometric  trans¬ 
formations.  There  is  a  proof  of  this  conjecture  for  some 
individual  cases  0. 

We  have  a  different  approach  because  the  proposed 
theory  generates  directed  acyclic  graphs  rather  than 
causal  sets.  These  graphs  have  partial  order  between 
pairs  of  adjacent  vertices,  so  causal  sets  appear  from 
them  via  transitivity  of  partial  order  on  directed  paths 
of  links.  In  addition  to  this,  the  vertices  of  the  graph 
have  a  measure  of  proximity  that  we  have  defined  as  the 
number  of  links  in  the  shortest  chain  of  links  regardless 
of  their  direction  l(v,w).  We  also  stick  to  the  paradigm 
of  “geometry  =  order  +  number” ,  but  we  have  to  define 
a  different  measure  of  volume  and  distance  because  the 
theory  does  not  have  the  required  properties  for  a  “faith¬ 
ful  embedding”  of  a  causal  set  into  a  spacetime  manifold. 
Here  is  a  review  on  the  CST  for  refferences  El- 

Spacetime  interval.  Now  it  is  time  to  introduce  an 
approach  to  define  a  covariant  measure  of  spacetime  in¬ 
terval  between  any  two  events.  We  start  with  a  definition 
for  a  spacelike  spacetime  interval. 
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FIG.  5.  A  schematic  example  for  calculation  of  spacelike 
spacetime  interval  SyW  between  events  u  and  w  by  finding 
direct  paths  l+(u,  v),  l+(v,  y),  l+(u,  w)  and  l+(w,  y),  that  min¬ 
imize  the  value  of  a  positive  spacelike  spacetime  interval. 


s2s(v,w)  =  min[(|4(u,uj)|  +  \l+(v,y)\)  x  (2a) 

u,y 

(I4(u,u)|  +  I4(w,t/)|)],  (2b) 

s.t.  u  £  J~(v)  D  J~(w),  (2c) 

y  £  J+  (v)  fl  J+  (w)  (2d) 

We  demonstrate  the  logic  behind  this  definition  with  a 
flat  spacetime  example  of  a  two-dimensional  square  lat¬ 
tice  shown  on  Figure  [5]  We  assume  that  all  directions 
of  edges  of  this  graph  go  from  bottom  to  top.  In  this 
example,  the  frame  of  reference  has  the  t-axis  along  one 
diagonal  of  the  square  lattice  and  the  x-axis  along  the 
other  diagonal  of  the  square  lattice.  The  coordinates  of 
all  events  are  the  following:  tu  =  0,  xu  =  0,  tv  =  2, 
xv  =  —2,  tw  =  4,  xw  =  4,  ty  =  6,  xy  =  2.  Then  the 
spacetime  interval  between  events  v  and  w  equals  to 


—  (Fw  &v)  if'W  tv)  (3a) 

-  (('T'U?  T  tyj)  ( Xy  T  ty ))  (3b) 

X  (( ty  Xy)  T  ( Xy]  tyj))  (3c) 

=  (2|4(u,  w) |  -  0)  x  (2|4(u,  y)\  +  0)  (3d) 

=  (|4(u,w)|  +  \l+(v,y)\)  (3e) 

x(|4(w,u)|  +  \l+(w,y)\)  (3f) 

=  8  *  4  =  32,  (3g) 


Definition  20  Spacelike  spacetime  interval  can  be  de¬ 
fined  for  any  two  spatially  separated  events  v  and  w  that 
do  not  have  any  directed  path  of  links  between  them  and 
have  a  joint  causal  future  J+{v)  fl  J+(w)  7^  0.  These 
events  always  have  a  common  past  J~{v)  fl  J~{w)  with 
at  least  the  initial  event  b.  A  positive  integer  value  of 
spacelike  spacetime  interval  equals  to: 


where  we  used  equations  |44,w)|  =  \l+(v,y)\  and 
144,  u)  |  =  \l+(w,y)\  for  a  case  of  flat  spacetime.  In  a 
general  case  of  non-flat  spacetime  we  simply  postulate 
Eq.  I3ell311  for  Eq.  [3a] 

Definition  21  Timelike  spacetime  interval  can  be  de¬ 
fined  for  any  two  timelike  separeated  events  v  and  w  that 


have  a  path  of  links  l+(v,w).  A  negative  integer  value  of 
timelike  spacetime  interval  equals  to: 


7 


St(v,  w)  =  -  max  [s2s(u,  y)]  ,  (4a) 

u,y 

s.t.  u,y  £  J+(v)  n  J~(w)  (4b) 

As  a  result,  we  have  defined  spacetime  interval  s2uw  for 
any  two  events  that  have  shared  future  purely  from  the 
structure  of  a  directed  acyclic  graph.  This  measure  is 
quantized,  and  it  is  a  multiplication  of  two  integers.  It 
is  either  positive,  negative  or  zero.  This  definition  for 
spacetinm  interval  is  different  from  that  in  the  causal  set 
theory  [TH  because  we  do  not  postulate  a  unit  volume 
for  each  event,  but  rather  provide  an  explicit  formula  for 
spacetime  interval. 

Time.  We  define  proper  time  on  a  world  line  or  any 
path  of  links  l+(v,w)  as  an  integral  sum  over  square 
roots  of  infinitesimal  timelike  spacetime  intervals  along 
the  world  line.  Since  the  spacetime  interval  between  ad¬ 
jacent  events  is  frequently  zero  while  it  is  non-zero  on  a 
larger  scale,  we  take  the  maximum  value  for  an  integral 
sum  over  all  segmentations  of  the  world  line  with  the 
length  of  each  segment  not  more  than  a  small  number  e, 
which  represents  the  size  of  a  clock  that  measures  time. 

Definition  22  Proper  time  along  a  path  l+(v,w)  is  de¬ 
fined  by  a  formula: 

At(1+(v,w))  =max  ^  \J~sly  (5) 

l+(u,y)€l+(v,w) 

Definition  23  Maximum  proper  time  between  any  two 
causally  related  events  with  multiple  possibilities  of  world 
lines  between  them  is  defined  by: 


Tmax(v,w)=  max  Ar(l+(v,w))  (6) 

l. |_  ( V,W ) 

Distance.  We  define  proper  length  or  rest  length 
along  a  chain  l(v,w )  between  causally  unrelated  events  v 
and  w  in  a  similar  way. 

Definition  24  Proper  length  between  two  causally  un¬ 
related  events  v  and  w  is  a  minimum  sum  over  square 
root  of  infinitesimal  spacelike  spacetime  intervals  along 
the  chain  l(v ,  w): 

A  L{l(v,w))  =  min  £  m 

l(u,y)€l(v,w ) 


Again,  this  definition  is  more  straightforward  than  the 
definition  in  the  causal  set  theory,  and  it  does  not  re¬ 
quire  special  consideration  for  3+  dimensional  spacetime, 
where  a  CST  definition  always  gives  the  value  of  two  as 
an  answer  0- 
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FIG.  6.  Dependence  of  the  radius  of  the  visible  Universe 
R-t  (u)  versus  the  maximum  proper  time  T-t  (v)  to  the  current 
epoch  for  an  analog  of  Model  A  in  four  dimensions  for  3128 
events  v  at  epoch  50  with  descendants  at  epoch  75.  This 
dependence  shows  the  edge  of  the  past  light  cone  that  is  linear 
for  small  times  T-t(v)  with  an  angle  consistent  with  the  speed 
of  light  c  =  1.  Every  point  is  an  average  over  all  events  v  for 
each  epoch  t. 


On  Figure  [6]  we  show  the  result  of  calculation  for  the 
edges  of  past  light  cones  in  four  dimensions.  As  expected, 
the  growth  of  the  average  radius  of  these  past  light  cones 
is  linear  to  the  maximum  proper  time  with  an  angle  ap¬ 
proximately  corresponding  to  the  speed  of  light  c  =  1. 
The  average  radius  also  decreases  to  zero  while  past  light 
cones  approach  the  initial  event  b. 

Inertial  frame  of  reference.  We  define  an  inertial 
frame  of  reference  in  the  following  way.  Let  us  consider  a 
region  of  spacetime  far  from  the  initial  event  b  and  away 
from  large-scale  curvature  anomalies.  Suppose  there  are 
two  events  v  and  w  with  the  shortest  path  between  them 
and  the  maximum  proper  time  Tmax(l* j_(u,w))  that  de¬ 
fines  a  world  line  for  the  origin  of  the  frame  of  refer¬ 
ence.  We  measure  time  by  the  proper  time  along  the 
path  Ijj _(v,w),  starting  from  the  event  v.  An  event  x  has 
the  same  time  as  an  event  y  €  l*L(v,  w)  if  the  event  y  max¬ 
imizes  spacelike  spacetime  interval  s2y  between  them. 
The  distance  between  two  events  with  equal  times  is  sim¬ 
ply  a  spatial  distance  between  them.  These  time  and 
distance  define  a  frame  of  reference  FRi.  We  call  it  an 
inertial  frame  of  reference.  In  general,  the  lines  of  equal 
time  and  distance  coordinates  cross  at  an  angle  that  de¬ 
pends  on  the  choice  of  the  frame  of  reference,  which  is 
set  by  a  pair  of  events  v  and  w.  We  keep  a  more  strict 
consideration  out  of  the  scope  of  this  paper. 

Stationary  points  in  FR\  are  spacetime  events  that 
have  world  lines  that  are  parallel  to  the  path 
i.e. ,  for  every  event  u  G  l*,(v,w)  the  maximal  spacelike 
spacetime  interval  to  the  events  on  these  parallel  world 
lines  is  the  same  up  to  a  small  value.  On  the  contrary,  a 


point  moving  with  a  constant  speed  in  the  frame  of  ref¬ 
erence  FRi  have  a  world  line  that  is  the  shortest  path 
between  two  events  on  world  lines  of  different  station¬ 
ary  points.  Suppose  a  moving  point  passes  by  an  event 
v,  so  its  world  line  is  along  a  path  l+(v,u).  The  angle 
between  l+(v,w)  and  l+(v,u)  defines  the  speed  of  the 
moving  point,  which  is  equal  to  the  spatial  distance  from 
u  to  l+(v,w)  divided  by  the  time  for  the  event  u.  See 
Figure  [7]  for  an  illustration. 

Doppler  effect.  The  initial  inspiration  for  this  the¬ 
ory  came  from  its  ability  to  demonstrate  mechanisms  for 
various  redshifts.  We  define  a  redshift  (blueshift)  as  an 
increase  (decrease)  of  maximum  proper  time  between  two 
points  at  the  edges  of  two  nested  light  cones. 

The  simplest  example  of  a  redshift  or  a  blueshift  is 
trivial  for  a  change  of  coordinates  from  one  inertial  frame 
of  reference  to  another,  taking  into  account  the  definition 
of  spacetime  interval.  This  is  a  demonstration  of  the 
Doppler  effect  [20j|. 

An  example  of  the  Doppler  effect  with  a  redshift  can 
be  demonstrated  within  the  considered  frame  of  reference 
FRi,  which  is  defined  by  events  v  and  w.  Suppose  there 
is  a  clock  that  moves  along  the  path  R(v,w)  and  the 
maximum  proper  time  between  v  and  w  is  tvw.  The 
clock  emits  information  at  every  tick.  Suppose  there  are 
k  equidistant  ticks  at  events  uq  £  l+{v,  w),  so  that  wo  =  v 
and  Wk  =  w.  The  frequency  of  emission  is  fvw  =  k/rvw. 
Every  tick  at  event  Wi  £  l+(v,w)  emits  a  future  light 
cone  J+(wi ),  so  that  the  points  at  the  edges  of  these 
nested  light  cones  move  away  from  the  origin  of  the  frame 
of  reference  FR±  with  a  maximal  speed  of  information 
propagation  or  the  speed  of  light.  If  we  take  any  event 
at  the  edge  of  the  last  light  cone  u  £  J+  (w)  then  we  can 
define  another  inertial  frame  of  reference  FR2  with  its 
origin  moving  along  the  shortest  path  u)  away  from 
the  origin  of  the  frame  of  reference  FR\.  The  maximum 
proper  time  between  v  and  u  is  tvu  .  Every  edge  of  a  light 
cone  J+(wi)  crosses  the  path  l+(v,  it),  so  the  tick  of  the 
clock  have  frequency  fvu  =  k/rvu  in  the  second  frame  of 
reference  FR2 .  The  redshift  happens  because  tvu  >  tvw  . 

The  proof  that  tvu  >  tvw  is  straigtforward.  By  defini¬ 
tion  of  spacetime  interval  for  a  flat  space:  s2w  =  —  Aa\b\, 
s2u  =  -4 a2b2,  and  |s2J  >  |s2J,  because  J+{u)  £ 
J+{w)  £  J+(v),  so  that  02  =  di  and  b2  >  b\.  Fi¬ 
nally,  from  the  definition  of  proper  time:  tvu  =  \J —s\u  > 
\/~suw  =  Tvw  ■  See  Figure  [7]  for  an  illustration. 

Relativistic  sum  of  velocities.  Derivation  of  the 
exact  Dopier  formula  and  a  relativistic  sum  of  two  ve¬ 
locities  can  be  done  in  a  similar  way  from  the  definition 
of  spacetime  interval,  proper  time  and  proper  distance. 
First,  we  define  velocity  versus  global  stationary  frame  of 
reference.  We  calculate  time  tvw  and  distance  xvw  coor¬ 
dinates  of  the  event  w  versus  the  event  v  in  this  frame  of 
reference  from  the  following  logic.  The  time  coordinate  of 
the  event  w  is  equal  to  two  terms:  tvw  =  tvy  +  tywi  where 

we  define  tvy  =  \j~s2y  =  yj 4 af  =  2oi  and  measure  tyw 
coordintate  along  a  light  like  spacitime  inteval  s2  =  0, 


t 


FIG.  7.  An  illustration  that  a  moving  object  w  that  is  em- 
miting  k  equidistant  light  cones  at  frequency  fvw  =  k/rvw 
is  seen  with  a  lower  frequency  of  emission  fvu  =  k/rVu  by  a 
faster  moving  object  u.  Events  y,  z  and  q  are  introduced  for 
calculation  of  coordinates  of  the  event  w. 


so  that  tyw  =  ^  y — syz  =  \j 4(6i  —  aq)2  =  b\  —  a\.  As 
a  result  tvw  =  &i  +  a±.  The  distance  coordinate  is  cal¬ 
culated  from  the  corresponding  space  like  spacetime  in¬ 
terval:  xvw  = 


=  i^/4(6i  -  a  1)2  =  bi  -  ai.  Now, 


'Vyw  —  Xyw  / tvw  • 

Finally,  we  get  a  genearal  formula  of  velocity  for  a 
world  line  with  a  timelike  spacetime  interval  s2  =  4 ab 
with  corresponding  paths  a  and  b  along  a  future  light 
cone  from  the  origin  of  the  frame  of  reference: 


b  —  a 
b  +  a 


(8) 


So  we  have  two  velocities  iq  =  (61  —  ai)/(b\  +  ai) 
and  v2  =  (b2  —  a2) / (b2  +  a2) ,  and  we  need  to  find 
a  formula  for  velocity  v'2  corresponding  to  velocity  v2 
in  the  frame  of  reference  FR\  that  is  moving  with  ve¬ 
locity  iq.  For  this  purpose,  we  need  to  calculate  co¬ 
ordinates  t'vu  and  x’vu  in  the  frame  of  reference  FR±. 
We  use  the  following  definition  of  spacetime  interval: 
s2vu  =  ~4a2b2  =  x2vu  -  t2vu  =  x,2u  -  t’2u.  And  we  use  the 
following  relation  between  the  coordinates  in  the  frame 
of  reference  FRy.  t'vu  =  t!vw  +t'wu  and  x'vu  =  x'wu.  Then 
we  use  the  dehnition  of  proper  time  t'vw  =  tvw  =  y/Aaibi 
and  the  definition  of  speed  of  light  x'wu  =  t'wu .  Finally, 
we  have  an  equation: 


-  4 a2b2  =  x,2u  -  (x/4aibi  +  x'wu)2 ,  (9) 

which  gives  a  solution  to  x'vu  =  ( a2b2  —  ai&i)/\/ai^i 
and  t'vu  =  (a2b2  +  a\b\) / y/aibi.  As  a  result,  we  calculate 
the  desired  velocity: 
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a2b2  —  aibi 
a2b2  +  a\bi 


(10) 


And  using  the  fact  that  a±  =  a2  we  come  to  a  well 
known  formula  of  a  relativistic  sum  of  velocities: 


1  —  v2  1  —  v'2  1  —  V\ 
1  +  V2  1  +  v'2  1  +  V\ 


(11) 


Finally,  we  get  the  exact  formula  for  the  Doppler  effect: 


f_  =  Tyw_=  jaih  _  I  l-v2 

f  rvu  V  a2b2  y  1  +  v'2  1  ’ 

Gravitational  redshift.  The  extremal  world  lines  at 
the  edges  of  light  cones  also  diverge  from  each  other  due 
to  the  local  curvature  of  spacetime  near  black  hole  sin¬ 
gularities.  In  order  to  see  this  effect  at  the  Planck  scale, 
we  show  an  example  of  two  future  light  cones  originat¬ 
ing  from  two  close  descendants  of  an  event  near  an  event 
horizon  of  a  black  hole,  see  Figure  [8]  The  edges  of  these 
two  future  light  cones  diverge  quickly  from  each  other 
as  it  happens  in  the  real  Universe  for  world  lines  of  two 
light  beams  originating  from  a  small  region  near  an  event 
horizon  of  a  black  hole. 

Cosmological  redshift.  There  is  global  curvature  in 
the  considered  theory  due  to  the  presence  of  the  initial 
event  b  at  the  center  of  the  graph.  This  global  curvature 
causes  continuous  divergence  of  the  edges  of  two  future 
light  cones  originating  from  two  close  events,  see  Figure [8] 
for  an  example.  The  same  logic  applies  to  the  past  light 
cones. 

Now  we  define  the  Hubble  constant  from  the  depen- 
dance  of  the  size  of  the  visible  Universe  and  the  time 
passage  from  that  epoch  M-  We  take  an  event  v  at 
epoch  n  and  all  events  in  its  past  light  cone  that  belongs 
to  epoch  n  —  t:  w  G  S-t(v)  =  J~(v)  H  E(n  —  t).  We  take 
the  maximum  proper  time  between  these  events: 


T-t(;v)=  max ■  r{l+(w,v)),  (13) 

weS-t(v) 

which  represents  the  amount  of  time  that  has  passed 
from  epoch  n  —  t  to  the  event  v  along  a  stationary  world 
line  originating  from  the  initial  event  b.  We  calculate  the 
maximum  diameter  between  the  events  from  the  space- 
time  region  S-t(v)  =  J~{y )  D  E(n  —  t)  as  a  maximum 
of  minimum  proper  distance  between  two  events  in  this 
region: 


D_t(u)  =  max  min  L(p(u,  y))  (14) 

u,y£S-t(v) 

We  now  define  a  radius  of  the  spacetime  region  S-t(v) 
as  a  maximum  distance  from  the  event  that  maximizes 


FIG.  8.  A  demonstration  of  strong  devergence  of  future  light 
cones  J+{x)  and  J+{y)  (light  blue  circles  and  yellow  circles) 
from  each  other  due  to  local  and  global  curvature  of  space- 
time,  regardless  the  fact  that  two  close  events  x  and  y  have 
a  common  preceeding  event  v.  Green  circles  show  an  epoch 
with  n  =  313  for  the  event  v.  Events  from  interiors  of  black 
holes  are  not  shown. 


proper  time  T_t(u).  At  a  large  scale  this  radius  is  a 
half  of  the  diameter:  R_t{v)  =  D_t(r;)/2.  In  flat  space 
R-t{v)  =  but  in  a  general  case  R-t{y )  /  T_t(u) 

becasue  of  spacetime  curvature.  Then  we  define  cosmo¬ 
logical  redshift  from  the  formula: 


1  +  z-t(v) 


T-t(v ) 
R-t(v) 


(15) 


The  redshift  z  is  zero  at  T  =  0,  and  it  is  positive  on 
a  large  scale  due  to  positive  overall  curvature  and  ex¬ 
pansion  of  space.  The  linear  term  in  the  Taylor  expan¬ 
sion  of  z(T)  is  related  to  the  Hubble  constant,  while  the 
quadratic  term  is  additionally  influenced  by  the  local  cur¬ 
vature:  <Z-t{v)  =  f(T-t{v))  ~  H(v)T-t(v),  where  H(v) 
is  an  approximate  value  of  the  Hubble  constant  around 
the  event  v.  It  is  measured  in  Planck  units.  Finally,  the 
value  of  the  Hubble  constant  at  epoch  n  is  an  average  for 
each  event  from  this  epoch:  H(n)  =  H(v). 

Calculation  of  the  Hubble  constant  from  the  ratio 
Z-t(v) /T-t(v)  at  t  ~  20  for  an  analog  of  Model  A  in  three 
dimensions  gives  an  approximate  value  of  H(n)  ~  — , 
which  is  in  line  with  the  current  estimation  for  H . 

Metric  tensor.  Now  we  introduce  an  embedding  of  a 
directed  acyclic  graph  Q  into  a  pseudo- Riemannian  man¬ 
ifold  M  with  a  mapping  r)  that: 

•  presersves  causality:  3 l+(v,w)  4$  r/(v )  -<  il(w) 


•  preserves  spacetime  interval  values  s2(v,w)  for  any 
events  v  and  w  with  common  future 


There  is  an  open  question  whether  this  mapping  77  ex¬ 
ists.  The  graph  Q  is  a  lattice,  so  there  is  no  problem 
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with  embedding  it  into  a  Euclidian  manifold.  The  crit¬ 
ical  question  is  whether  a  pseudo-Riemannian  signature 
metric  tensor  g exists  that  is  consistent  with  all  causal 
relations  and  spacetime  intervals. 

The  mapping  ij  is  also  non-unique  at  least  for  small 
deformations  of  A i  at  the  scale  smaller  than  a  link  of 
the  graph.  There  is  a  proven  theorem  that  a  Lorentzian 
geometry  A4  follow  from  causal  order  only  up  to  a  lo¬ 
cal  conformal  factor  [22|,  [23[,  [24[,  j5H[.  Unit  volume 
y/— gdx 4  fixes  this  local  conformal  factor  in  the  main¬ 
stream  causal  set  theory.  The  CDFT,  on  the  contrary, 
fixes  timelike  and  spacelike  spacetime  intervals  between 
the  vertices  of  a  directed  acyclic  graph. 

In  order  to  move  forward,  we  formulate  a  variant  of 
the  central  conjecture  for  the  causal  discrete  field  theory: 

Every  directed  acyclic  graph  Q  from  the 
causal  discrete  field  theory  can  be  uniquely 
embedded  into  a  pseudo-Riemannian  man¬ 
ifold  up  to  an  isometric  transformation: 

Mi(g)^M2(G) 

We  propose  an  approximate  way  to  calculate  the  met¬ 
ric  tensor  gpl,(v)  for  each  event  v  from  the  geometry  of 
its  vicinity  V ( v ,  e)  that  is  e  links  apart  by  minimizing  the 
fitting  error  of  the  local  metric  tensor  to  reproduce  local 
spacetime  intervals. 

First,  we  assume  some  vector  coordinates  for  an  n- 
dimensional  case  with  the  center  at  the  initial  vertex  b. 
We  use  unit  length  basis  vectors  along  n  chosen  graph 
links  originating  from  the  initial  vertex  b.  Any  choice 
of  n  links  out  of  n  +  1  is  possible  because  the  graph  is 
self-similar.  The  coordinates  of  any  event  appear  from 
parallel  translations  of  this  unit  length  basis  vectors. 
This  definition  of  coordinates  is  proper  for  simple  self¬ 
similar  graph  geometries  like  cubic  lattice  or  generalized 
diamond  lattice.  Some  exotic  self-similar  graphs  may 
need  a  different  approach  with  non-unique  coordinates 
for  the  same  events,  e.g.,  for  a  finite  graph  like  a  fullerene 
molecule  or  a  graph  with  some  compactified  finite  dimen¬ 
sions  like  a  nanotube  molecule. 

Second,  we  take  all  events  from  the  vicinity  of  the 
givent  event:  w  £  V(v,e),  and  calculate  coordinate  dif¬ 
ferentials  between  the  events  dx%w  =  —  x%  and  actual 

spacetime  intervals  between  them  s^w . 

Finally,  we  minimize  square  errors  between  a  bi-linear 
combination  of  coordinate  differentials  and  actual  space- 
time  intervals: 

slw  =  SwK  e)dx^wdxuvw  +  evw,  (16a) 

s.t.  e„w  —>  min  (16b) 

wE.V(v,e) 

As  a  result,  we  calculate  an  e-averaged  metric  tensor 
g^v,  where  e  is  the  size  of  the  averaging  region  measured 
in  the  number  of  links  from  the  chosen  event  v. 

Spacetime  curvature.  We  propose  to  do  an  approx¬ 
imate  calculation  of  all  relevant  elements  of  Riemannian 


geometry.  We  can  take  approximate  values  of  metric  ten¬ 
sor  at  each  event  g^(v,  e)  and  calculate  a  twice  differen¬ 
tiable  interpolation  of  it  beyond  the  vertices  of  the  graph. 
It  allows  to  calculate  partial  derivatives  of  the  metric  ten¬ 
sor: 


d  8 2 

dxk9ij~9ij,k  and  gxkgxl9ij  ~  9ij,kl 


(17) 


And  then  combine  them  into  the  Christoffel  symbols: 


1  ( _d_  _d_  d  \ 

kij~  2  \d^9H  +  dtf9kj  ~ 

=  l(9kij+9ki,i-9H,k) 


(18a) 

(18b) 


T™  =  gmkTkij  =  X-gmk (gkiJ  +  gkjti  -  gijtk)  (19) 
Finally,  the  Riemannian  curvature  tensor  equals  to 


r)  f) 

jdI  _  _  -pi  _  _  W  |  -pi  -ps 

11  ijk  —  gxjL  ik  dxkL  b'  "t_  1  JS1  ik 
Lowering  the  indexes  Rsijk  =  gsiRl gk 


_ Y'l  ps 

1  ks ij 

we  get 


(20) 


Rsijk  2  (9sk,ij  T  9ij.sk:  fj-s  j/ik  fjik.sj )  T  (21a) 

+gnp(TnijTpsk  -  T"ifcrW)  (21b) 

Finally,  Ricci  curvature  tensor  and  the  Ricci  curvature 
equal  to 


Rij  =  glmRiijm  and  JZ  =  g^RtJ  (22) 

We  have  described  a  way  to  calculate  approximate  met¬ 
ric  tensor  and  approximate  curvature  tensor  appealing  to 
continuous  Riemannian  geometry.  It  is  helpful,  for  ex¬ 
ample,  to  calculate  large-scale  curvature  of  the  resulting 
spacetime.  However,  it  is  still  interesting  to  find  a  way  to 
define  discrete  analogs  of  these  objects  purely  from  the 
structure  of  a  directed  acyclic  graph. 

Energy  and  momentum.  Now  we  introduce  a  trivial 
definition  of  energy  and  momentum.  There  is  already  a 
metric  tensor  gjj  in  this  theory.  This  tensor  field  defines 
all  objects  within  Riemannian  geometry  on  a  manifold, 
including  the  Einstein  tensor: 


Gij  =  Rij  -  -Rg^  (23) 

The  Einstein  tensor  is  symmetric  Gij  =  Gji  and  it 
is  divergenceless  =  Gl.j  =  0.  The  stress-energy 

tensor  of  matter  and  the  metric  tensor  have  the  same 
property,  so  we  can  postulate  the  stress-energy  tensor  to 
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FIG.  9.  Dependence  of  the  size  of  the  whole  Universe  L-t(v) 
versus  the  maximum  proper  time  T-t(v)  to  the  current  epoch 
for  an  analog  of  Model  A  in  three  dimensions  for  1747  events 
v  at  epoch  200  with  descendants  at  epoch  250.  Every  point 
is  an  average  over  all  events  v  for  each  epoch  t. 


FIG.  10.  Dependence  of  the  size  of  the  whole  Universe  L-t(v) 
versus  the  maximum  proper  time  T-t(v)  to  the  current  epoch 
for  an  analog  of  Model  A  in  four  dimensions  for  3128  events 
v  at  epoch  50  with  descendants  at  epoch  75.  Every  point  is 
an  average  over  all  events  v  for  each  epoch  t. 


be  equal  to  the  Einstein  tensor  up  to  a  cosmological  term 
given  the  Planck  units  assumption  of  G  =  1  and  c  =  1: 

Tij  =  —  (Gij  +  Agij )  (24) 


for  an  analog  of  Model  A  in  three  and  four  dimensions 
because  the  two-dimensional  case  is  degenerate.  We  es¬ 
timate  the  size  of  the  whole  Universe  at  each  epoch  n  —  t 
from  the  radius  of  the  visible  Universe  R-t(v)  and  its 
angular  size  </>_t(u)  in  the  embedding  space  of  the  graph: 


A  more  strict  approach  requires  an  understanding  of 
the  local  dynamics  of  discrete  fields  on  the  graph.  We  use 
the  Einstein  field  equations  to  reconstruct  the  energy- 
momentum  tensor  because  gravitational  field  emerges 
selfconsistently  without  equations  of  matter.  We  develop 
the  theory  from  pure  geometry  of  the  graph  and  the  prop¬ 
agation  rule. 

The  role  of  the  cosmological  term  A  is  still  not  clear. 
This  term  can  also  be  variable,  which  leads  to  the  trans¬ 
fer  of  energy  and  momentum  between  the  cosmological 
term  and  matter.  The  implications  of  a  variable  cos¬ 
mological  term  for  evolution  of  the  scale  factor  a(t)  was 
studied  in  literature  with  different  assumptions  for  decay 
laws  {26).  Our  favorable  hypothesis  is  to  use  it  to  offset 
the  global  curvature  of  spacetime.  The  idea  is  to  have 
an  energy-momentum  tensor  T. \j  that  is  only  due  to  the 
local  dynamics  of  fields  on  the  graph.  As  a  result,  the 
cosmological  term  would  be  positive  and  proportional  to 
the  inverse  square  of  epoch  number  n:  A  =  7 Z(n)  ~  , 

which  is  in  line  with  the  current  estimation  for  A. 

Inflation.  Computable  Model  A  does  not  reproduce 
the  inflation  scenario  right  after  the  Big  Bang.  It  is  an 
open  question  to  check  this  matter  for  non-computable 
models.  The  inflationary  expansion  of  the  Universe  {27] 
was  proposed  to  explain  the  uniformity  of  cosmic  back¬ 
ground  radiation  [281]  that  has  not  interacted  with  matter 
since  the  moment  of  the  photon  decoupling  event  [29| . 

We  calculate  all  spacetime  intervals  and  proper  times 


L-tiy)  =  ——R_t(v)  (25) 

<p-t(y) 

The  size  of  the  Universe  L_t(u)  depends  linearly  from 
the  proper  time  T_t(v),  see  Figures  IT71  and  [TUI  In  four 
dimensions:  T  sa  —0.391  and  L  ss  27t  x  2.1(n  —  t). 

We  also  compare  our  numeric  calculations  with  a  con¬ 
tinuous  limit.  In  this  case  the  edge  of  the  past  light  cone 
is  described  by  an  equation  in  polar  coordinates  of  the 
angle  <j)  and  the  overall  epoch  number  t: 


, ,  cadt 
d<p  = - 


(26) 


where  we  have  introduced  a  rate  of  time  flow  which  cor¬ 
responds  to  an  angle  of  the  edge  of  the  past  light  cone  in 
the  embedding  space  of  the  underlying  graph:  a  =  jjj . 
With  a  constant  speed  of  light  c  =  1  and  a  constant  rate 
of  time  flow  a,  the  angle  coordinate  for  the  edge  of  the 
past  light  cone  starging  at  epoch  t0  will  depend  on  the 
past  epoch  t\  in  the  following  way: 


Atj)  =  a  ln(— )  (27) 

t  i 

This  formula  shows  that  if  one  looks  back  for  long 
enough,  then  the  whole  Universe  can  be  seen  an  infinite 
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number  of  times  in  the  continuous  limit.  If  in  practice 
the  backward  propagation  stays  within  a  limited  angle 
because  the  rate  of  time  flow  a  is  much  lower  for  smaller 
epochs,  then  a  distant  observer  at  the  current  epoch  to 
will  perceive  the  Universe  expanding  exponentially  fast 
close  to  the  Big  Bang  event. 

We  show  our  calculations  for  the  angular  size  of  the 
visible  Universe  on  Figures  |TT]  and  [12]  We  make  the  best 
fit  for  03 d  ~  xjg  in  three  dimensions  and  a^d  ~  <yg  in 
four  dimensions,  which  represents  the  angle  for  the  edges 
of  past  light  cones  at  large  epochs  and  approximately 
equals  to  a  ~  2g£  x  It  would  be  interesting  to  find 
the  exact  value  of  a  for  Model  A  with  large  n. 

The  simulated  dependence  in  three  dimensions  is  close 
to  the  continuous  limit  with  slightly  slower  expansion 
during  the  early  epochs  and  marginally  faster  expansion 
at  the  middle  of  the  plot.  The  simulated  dependence  in 
four  dimension  fits  the  line  of  the  continuous  limit  even 
better,  but  it  has  only  50  calculated  epochs.  This  com¬ 
putational  result  suggests  that  Model  A  generates  a  sim¬ 
ple  cosmological  model  [l^l  with  a  variable  cosmological 
term  A  ~  ~  -^y,  which  avoids  the  problem  of  small 

fundamental  constants  [3l|,  [H}.  The  resulting  metrics 
depends  on  the  epoch  t  similar  to  the  Milne  model 
but  with  positive  curvature: 


FIG.  11.  Dependence  of  the  angular  size  of  the  visible  Uni¬ 
verse  (p~t(v )  from  epochs  t:  simulated  points  (circles)  and 
continuous  limit  with  a^d.  ~  y^-g  (solid  line).  The  calculation 
is  done  for  an  analog  of  Model  A  in  three  dimensions  for  1747 
events  v  at  epoch  200  with  descendants  at  epoch  250.  Every 
point  is  an  average  over  all  events  v  for  each  epoch  t. 


ds 2  =  — c2dt 2  +  —z  (  ^  +  r2dff2^  (28) 

cr  \  1  —  rz  ) 

The  inflationary  scenario  requires  that  simulated  an¬ 
gles  for  early  epochs  are  located  much  lower  than  the 
theoretical  line  for  the  continuous  limit.  An  alternative 
explanation  for  a  homogeneous  Universe  at  the  photon 
decoupling  event  is  that  the  Universe  was  causally  con¬ 
nected  at  early  epochs.  In  the  continuous  limit  of  four 
dimensions,  the  edges  of  the  past  light  cone  can  do  three 
circles  around  the  Big  Bang  event  during  to  ~  8  x  106° 
Planck  units  of  time,  and  the  visible  Universe  coincided 
with  the  whole  Universe  at  the  moment  of  A =  tt  at  the 
age  of  the  Universe  of  ii  ~  7.3  years. 

CPT  invariance.  Symmetry  plays  an  essential  role 
in  theoretical  physics.  So  far,  we  have  discussed  trans¬ 
lational  symmetry,  which  follows  from  the  fact  that  each 
vertex  has  the  same  local  dynamics  due  to  a  propagation 
rule.  We  also  conjectured  that  rotational  symmetry  and 
isotropy  of  space  follow  from  isotropic  propagation  rules 
for  some  types  of  graphs,  e.g.,  for  generalized  diamond 
crystal  lattices.  Now  we  turn  to  another  fundamental 
symmetry  for  matter  and  anti-matter,  which  is  described 
by  CPT  invariance. 

Time  reversal  symmetry  changes  the  direction  of  all  di¬ 
rected  links  on  the  graph:  ipij  =  +c  4=>  ipij  =  —c.  Parity 
symmetry  makes  a  mirror-image  of  the  graph.  Charge 
symmetry  permutes  integer  values  on  edges:  i'  =  C(i), 
where  C(i)  is  an  integer  defined  function  that  exchanges 
some  pairs  of  integers  or  keeps  some  integers  unchanged. 

Each  of  these  symmetry  transformations  matches  the 
graph  with  itself  and  changes  the  pattern  of  discrete  fields 


Epoch  (t) 


FIG.  12.  Dependence  of  the  angular  size  of  the  visible  Uni¬ 
verse  )  from  epochs  t:  simulated  points  (circles)  and 

continuous  limit  with  aid  ~  gy.  (solid  line).  The  calculation 
is  done  for  an  analog  of  Model  A  in  four  dimensions  for  3128 
events  v  at  epoch  50  with  descendants  at  epoch  75.  Every 
point  is  an  average  over  all  events  v  for  each  epoch  t. 


on  it.  We  call  two  propagation  rules  symmetrical  if  and 
only  if  for  each  pattern  on  a  graph  generated  by  one  prop¬ 
agation  rule,  there  is  a  symmetric  pattern  on  a  graph  gen¬ 
erated  by  the  other  propagation  rule.  C,  P  or  T  symmet¬ 
ric  propagation  rules  can  appear  from  the  initial  prop¬ 
agation  rule  by  applying  C(i)  to  integers  on  incoming 
and  outgoing  edges,  applying  a  non-parity  permutation 
of  edges  or  by  changing  the  direction  of  all  incoming  and 


13 


outgoing  edges. 

Each  of  the  three  symmetries  allows  for  a  non-trivial 
change  of  a  propagation  rule  so  that  the  new  rule  is  sub¬ 
stantially  different  from  the  initial  one.  A  notable  case 
is  when  the  application  of  all  three  symmetries  simulta¬ 
neously  creates  the  same  propagation  rule.  We  call  this 
type  of  propagation  rules  as  CPT  invariant  propagation 
rules. 

We  call  a  propagation  rule  violating  symmetry  if  and 
only  if  the  application  of  this  symmetry  to  all  its  patterns 
on  graphs  substantially  changes  the  resulting  propaga¬ 
tion  rule. 

Model  A  and  Model  A-D  on  Figure  |T|  are  examples 
with  CPT  invariant  non-trivial  propagation  rules.  Model 
A-D  also  violates  C,  P,  and  T  symmetries  separately,  be¬ 
cause  all  four  vertices  in  sector  B  do  not  have  symmetrical 
counterparts  for  all  three  simple  symmetries. 

The  most  exciting  types  of  propagation  rules  are  non¬ 
trivial  non-computable  CPT  invariant  isotropic  propaga¬ 
tion  rules  with  violation  of  separate  C,  P  or  T  symmetries 
because  they  are  most  similar  to  what  happens  in  the  real 
Universe.  However,  these  propagation  rules  are  hard  to 
compute,  which  complicates  the  analysis. 

An  interesting  topic  is  the  CPT  theorem  stating  that 
any  Lorentz  invariant  local  quantum  field  theory  with  a 
Hermitian  Hamiltonian  must  have  a  CPT  symmetry  [34| . 
If  there  is  an  analog  of  this  theorem  for  the  proposed 
theory,  then  CPT  invariant  propagation  rules  must  be 
the  key  for  finding  a  quantum  mechanical  description  for 
CDFT  models. 

Elementary  particles.  There  are  still  no  particles  in 
the  overall  picture.  The  proposed  theory  is  self-consistent 
so  that  it  does  not  require  special  consideration  of  matter 
or  elementary  particles.  However,  propagation  rules  may 
have  one  more  symmetry  of  charge  conservation,  which 
together  with  CPT  invariance  may  be  important  for  the 
introduction  of  elementary  particles  in  this  theory. 

Model  A-D  on  Figured] is  a  simple  example  of  charge 
conservation:  graph  edges  with  integer  2  are  currents 
with  positive  charge  and  graph  edges  with  integer  3  are 
currents  with  negative  charge,  and  the  propagation  rule 
allows  only  uninterrupted  flow  of  these  charges  in  addi¬ 
tion  to  creation  and  annihilation  of  pairs  of  currents  with 
opposite  charges. 

Some  propagation  rules  can  also  have  conservation  of 
multiple  types  of  charges,  for  example,  by  having  mul¬ 
tiple  pairs  of  currents  like  in  Model  A-D.  Also,  partial 
charge  conservation  is  possible,  like  in  the  real  Universe. 
For  example,  there  can  be  strict  conservation  of  “elec¬ 
tric”  charge  when  the  number  of  odd/even  integer  value 
incoming  edges  equal  to  the  number  of  odd/even  integer 
value  outgoing  edges  for  each  vertex  except  for  “neutral” 
currents  like  integer  1  in  Model  A-D.  At  the  same  time, 
there  can  be  pairs  of  odd/even  integer  values  j  <t=>  C(i) 
that  are  are  “mixed”  without  parity  violation  for  some 
rear  types  of  vertices. 

Mass  and  inertia.  We  now  suggest  a  mechanism  of 
how  masses  might  appear  in  CDFT  models.  We  already 


postulated  that  every  field  configuration  can  happen  with 
a  likelihood  that  is  defined  by  a  constraint  functional 
C-p in  Eq.  [T]  In  quantum  physics,  the  principle  of 
least  action  defines  which  field  configurations  or  paths  of 
particles  are  more  likely  than  others.  We  also  follow  that 
approach. 

The  likelihood  of  taking  an  infinite  size  pattern  from 
an  infinite  uncountable  set  is  poorly  defined,  so  we  start 
from  the  likelihood  of  local  configurations  for  valid  types 
of  vertices.  We  postulate  the  following  “principle  of  least 
action” : 

Sequential  growth  of  vertices  is  equally  likely 
in  terms  of  the  choice  of  all  possible  locations 
and  types  for  the  next  new  vertex. 

This  principle  makes  it  easy  to  construct  patterns  on 
self-similar  graphs  with  computable  propagation  rules 
by  randomly  choosing  which  vertex  of  which  valid  type 
to  add  next  to  the  growing  set  of  already  enumerated 
vertices.  Non-computable  propagation  rules  addition¬ 
ally  need  to  have  a  non-computable  “oracle”  that  knows 
which  types  of  vertices  are  valid  locally  given  global  con¬ 
straints  in  the  future. 

The  physical  meaning  of  this  principle  is  the  following. 
A  local  observer  does  not  care  about  the  global  likeli¬ 
hood  of  future  patterns.  All  that  she  cares  is  that  there 
are  some  possible  splits  of  her  Universe  on  the  next  ver¬ 
tex  of  her  world  line,  and  she  can  find  herself  in  any  one 
of  them  with  equal  likelihood.  This  postselection  prob¬ 
ability  calculation  resembles  the  famous  sleeping  beauty 
paradox  when  an  observer  assigns  equal  subjective  prob¬ 
abilities  for  indistinguishable  future  outcomes. 

The  proposed  postselection  of  Universe  splits  produces 
a  tremendous  difference  in  the  likelihood  of  valid  patterns 
like  it  is  for  possible  virtual  paths  of  elementary  parti¬ 
cles  in  quantum  field  theories.  Vertices  with  “charged” 
edges  are  expected  to  have  a  lower  likelihood  on  average 
because  charge  conservation  is  an  additional  constraint 
that  needs  to  be  satisfied.  So  sequential  growth  might  be 
slower  statistically  near  conserving  charges,  which  should 
lead  to  the  emergence  of  persistent  craters  with  positive 
Riemannian  curvature  around  them.  That  is  how  parti¬ 
cle  and  antiparticle  charges  might  acquire  equal  positive 
masses.  Gravity  might  be  a  result  of  a  massive  particle 
being  placed  in  curved  spacetime  around  another  massive 
particle. 

Another  direction  to  explore  is  whether  the  sequential 
growth  of  vertices  allows  for  charges  moving  with  a  con¬ 
stant  speed  versus  the  global  frame  of  reference.  It  is  nat¬ 
ural  to  expect  stable  symmetric  craters  around  charged 
world  lines  going  perpendicular  to  the  surface  of  sequen¬ 
tial  growth.  The  question  is  whether  locally  equiproba- 
ble  sequential  growth  can  generate  stable  skewed  craters 
around  charged  world  lines  going  at  an  angle  to  the  sur¬ 
face  of  sequential  growth.  This  would  mean  that  free 
charges  have  not  only  mass  but  also  inertia. 

Both  ideas  for  mass  and  inertia  of  free  moving  charges 
are  entirely  speculative  and  up  to  future  confirmation  be- 
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cause  charge  conservation  implies  non-computable  prop¬ 
agation  rules  that  are  hard  to  model  computationally. 
The  theory  allows  for  a  huge  number  of  field  configura¬ 
tions  with  the  possibility  to  reproduce  any  complex  inter¬ 
action  of  particles  and  spacetime  metrics  around  them,  so 
the  question  is  whether  physically  meaningful  field  con¬ 
figurations  are  more  likely  due  to  equiprobable  sequential 
propagation. 

Entropy.  Locally  CPT  invariant  patterns  on  a  graph 
still  have  an  essential  global  asymmetry  versus  time  re¬ 
versal,  because  we  require  only  one  initial  vertex  and 
equal  likelihood  for  local  types  of  vertices  during  sequen¬ 
tial  growth.  On  the  contrary,  there  are  a  lot  of  terminal 
vertices,  which  are  singularity  points  of  black  holes.  So 
time  reversal  of  this  kind  of  Universe  will  lead  to  only  one 
giant  black  hole  and  an  infinite  number  of  “white  holes” , 
which  are  CPT  analogs  of  black  holes.  In  the  real  Uni¬ 
verse,  we  do  have  signs  of  existence  for  black  holes  [35|. 
while  “white  holes”  remain  only  a  theoretical  object. 

This  global  time  asymmetry  may  explain  the  arrow  of 
time  and  the  second  law  of  thermodynamics  that  forbids 
the  decrease  of  entropy  in  a  closed  system.  In  the  pro¬ 
posed  theory  the  time-reversed  Universe  indeed  collapses 
to  a  single  state  of  the  terminal  point  of  a  giant  black  hole, 
so  the  entropy  of  the  whole  Universe  eventually  drops  to 
the  least  possible  value  of  zero  at  the  Big  Bang  event. 

Here  is  an  example  of  a  thought  experiment  to  calcu¬ 
late  entropy  within  the  proposed  theory.  We  consider 
a  case  of  a  “stationary”  Universe  in  two  dimensions  for 
simplicity.  We  use  a  nanotube  molecule  type  of  graph 
with  one  class  of  the  edges  of  the  graph  going  along  the 
nanotube.  This  graph  has  one  infinite  dimension  along 
the  nanotube  and  one  confined  dimension  perpendicu¬ 
lar  to  it.  The  diameter  of  the  nanotube  should  be  large 
enough  to  dismiss  the  effects  of  circular  boundary  con¬ 
ditions.  Suppose  we  have  n  vertices  in  each  circle  of  the 
nanotube  that  go  with  a  period  of  four  given  the  number 
of  different  vertices  of  a  hexagon  along  the  nanotube. 

We  propose  to  seed  initial  vertices  on  the  graph  as  a 
circle  chain  of  valid  types  of  vertices  around  the  nan¬ 
otube.  Then  we  propose  to  run  sequential  growth  in  one 
of  the  directions  along  the  nanotube  from  the  initial  ver¬ 
tices  for  long  enough  to  get  statistical  equilibrium.  The 
direction  of  propagation  will  correspond  to  a  global  time 
arrow  in  this  stationary  Universe. 

We  propose  to  use  a  CPT  invariant  charge  conserving 
propagation  rule  so  that  the  flow  of  charge  through  the 
nanotube  is  constant  given  the  convention  that  flow  of 
charge  c  in  the  reverse  time  direction  corresponds  to  a 
flow  of  charge  C(c)  in  the  global  time  direction. 

Definition  25  A  microstate  |<fq)  is  a  set  of  open-ended 
graph  edges  during  sequential  growth  without  standalone 
cavities  that  eventually  collapse  into  black  holes. 

Each  step  of  sequential  growth  adds  one  vertex  by  con¬ 
necting  it  to  one  or  more  existing  open-ended  graph  edges 
and  forming  zero  or  more  new  open-ended  edges  from 
the  new  vertex.  We  exclude  cavities  from  consideration 


because  they  are  not  causally  related  to  the  remaining 
part  of  the  graph  any  more,  and  therefore  they  should  be 
treated  as  separate  microstates. 

By  running  sequential  growth,  we  may  reach  a  statis¬ 
tical  equilibrium  and  calculate  the  probability  for  each 
microstate  |<Iq).  The  number  of  all  possible  microstates 
is  infinite  because  there  is  an  infinite  number  of  elongated 
ellipses  around  the  nanotube.  However,  we  expect  that 
circular  shape  microstates  should  dominate  the  probabil¬ 
ities  so  that  infinite  sums  over  all  microstates  will  con¬ 
verge. 

Definition  26  A  macrostate  |  T)  is  a  set  of  microstates 
|<Eq)  that  have  the  same  property  or  symmetry. 

We  consider  two  variants  of  macrostates  in  our  exam¬ 
ple.  First  macrostate  has  a  zero  flow  of  charge:  | \Ho)  - 
Second  macrostate  has  a  unit  flow  of  charge:  I'Ll).  Since 
we  assume  a  charge-preserving  propagation  rule,  all  mi¬ 
crostates  will  stay  within  their  corresponding  macrostate 
during  propagation. 

Now  we  can  calculate  entropy  for  each  macrostate  with 
charge  q  in  equilibrium  with  a  Planck  value  for  the  Bolz- 
mann  constant  kp  =  1: 

^l^9>  =  —  5Z  Pi^Pi  (29) 

1 4-del**? 

We  anticipate  that  the  marcostate  with  charge  | has 
higher  entropy  than  the  macrostate  without  charge  |  H/o) , 
because  due  to  charge  conservation  constraints  fewer  mi¬ 
crostates  |<P;)  will  have  high  enough  propabilities  pi. 

The  delta  between  the  two  entropies  corresponds  to 
the  entropy  of  a  massive  particle  versus  a  vacuum  state: 

AS1  =  S,|4,1)  —  S|*0>  >  0  (30) 

This  entropy  A S  adds  another  way  in  addition  to  the 
Einstein  field  equations  in  Eq.  [24]  to  estimate  mass  by 
analogy  with  the  Bckenstein-Hawking  entropy  of  a  sym¬ 
metric  Schwarzschild  black  hole,  which  is  proportional  to 
the  square  of  its  mass  in  four-dimensional  spacetime: 

SBh=4ttM2,  (31) 

where  M  is  the  mass  of  the  black  hole  in  Planck  units  of 
mass  mp  =  ^  =  1. 

It  is  an  interesting  open  question  to  prove  the  power 
law  dependence  of  entropy  in  Eq.  [30]  from  geometrically 
defined  mass  in  Eq.  [24]  in  various  dimensions  and  to  find 
the  exact  coefficient  in  the  equation  similar  to  Eq.  1311 

Quantum  field  theory.  The  most  challenging  task 
for  the  CDFT  is  to  explain  quantum  effects  within  a 
causal  and  discrete  paradigm.  One  alternative  is  to  de¬ 
fine  a  continuous  field  theory  that  has  a  field  9  on  some 
pseudo-Riemannian  n-dimensional  manifold  and  an  ac¬ 
tion  term  ,  so  that  this  field  theory  allows  for  a  quan¬ 
tization  constraint  that  makes  a  correspondence  between 
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the  fields  [0]  <t=>  [ip\  and  the  action  terms  elS^  +>  C-p,b [V>] 
for  a  certain  self-similar  graph. 

So  far,  we  have  only  been  able  to  outline  a  way  to  find 
a  correspondence  between  a  discrete  field  of  directions 
of  links  [^]  and  a  continuous  tensor  field  for  metrics  g 
without  even  construction  of  the  action  for  the  gravity 
part  of  the  theory. 

We  propose  to  search  for  a  variant  of  a  CDFT  model 
that  allows  for  a  continuous  quantum  field  approximation 
that  is  consistent  with  the  Standard  Model  of  elementary 
particles  at  large  scale  and  low-energies  in  flat  spacetime 
far  away  from  the  Big  Bang  event. 

Some  progress  has  been  made  recently  with  another 
approach  that  shows  how  a  deterministic  cellular  au¬ 
tomata  model  can  demonstrate  quantum  field  theory 
phenomena.  Gerard  ‘t  Hooft  has  studied  a  simple  cellular 
automata  model  that  is  equivalent  to  string  theory  in  1+1 
dimensions  [36[.  However,  a  cellular  automata  model  de¬ 
fines  propagation  of  information  on  separate  slices  of  time 
and  therefore  has  a  problem  describing  covariant  objects 
in  curved  spacetime. 

As  the  next  step  for  future  research,  we  propose  to 
develop  these  ideas  further  adding  the  many-worlds  in¬ 
terpretation  as  an  essential  ingredient,  where  a  local  pat¬ 
tern  around  one  event  unitarily  evolves  to  a  mix  of  local 
patterns  around  another  event. 

Many-worlds  interpretation.  The  CDFT  is  con¬ 
sistent  with  the  many-worlds  interpretation  of  quantum 
mechanics  that  took  its  roots  from  Erwin  Schrodinger’s 
ideas  [37:|,  was  first  described  by  Huge  Everett  [38l.  (doj 
and  was  later  popularized  by  Bryce  DeWitt  [4f|.  The 
Multiverse  in  the  CDFT  contains  all  copies  of  Universes, 
each  pair  of  which  has  some  maximal  common  area 
around  the  Big  Bang  event.  The  Multiverse  splits  any 
time  when  the  propagation  rule  allows  for  more  than  one 
configuration  of  integer  fields  and  edge  directions  for  the 
same  variant  of  the  past.  This  ambiguity  can  happen 
either  due  to  multiple  valid  types  for  the  next  vertex  in 
sequential  growth  or  due  to  the  application  of  a  propaga¬ 
tion  rule  in  a  different  order  for  this  vertex  and  its  direct 
neighbors. 

The  CDFT  has  the  following  variant  of  the  many- 
worlds  interpretation  of  quantum  mechanics: 

•  Infinite  uncountable  number  of  Universes 

•  Finite  number  of  Universes  of  a  given  size 

•  Splitting  happens  only  in  causal  future 

•  No  interaction  between  different  Universes 

•  The  whole  Multiverse  is  generated  by  a  determin¬ 
istic  propagation  rule 

•  The  Multiverse  is  computable  for  computable  prop¬ 
agation  rules 

•  The  Multiverse  may  be  hard  to  compute  for  non- 
computable  propagation  rules 


The  resulting  Multiverse  does  not  only  have  multi¬ 
ple  variants  of  Universes,  but  it  also  allows  for  all  four 
non-exclusive  possibilities  that  are  considered  in  litera¬ 
ture  [H}.  First,  sub-universes  may  be  different  regions 
in  space.  Indeed,  some  events  may  not  have  a  shared 
future,  and  they  may  have  a  shared  past  only  very  close 
to  the  Big  Bang  event.  Second,  fundamental  properties 
may  change  with  time.  We  postulated  that  the  speed 
of  the  fastest  information  propagation  is  constant,  and 
we  derived  time  and  spatial  distance  from  that.  How¬ 
ever,  some  other  physical  constants  may  change  slowly, 
for  example,  the  cosmological  term  A  and  the  Hubble 
constant  H  depend  on  the  global  curvature  of  spacetime. 
Third,  the  Universe  can  have  spacetime  regions  with 
different  properties.  Information  patterns  on  a  graph 
may  undergo  phase  transitions  while  the  global  curva¬ 
ture  changes  so  that  the  propagation  of  information  pat¬ 
terns  on  the  graph  happens  in  different  ways.  Finally, 
the  real  Universe  can  be  a  subspace  of  a  higher  dimen¬ 
sional  CDFT  model  with  some  dimensionality  reduction 
to  a  four-dimensional  spacetime. 

Apart  from  quantum  mechanics,  the  discrete  version 
of  the  many-worlds  interpretation  provides  an  interesting 
view  on  the  von  Neumann-Morgenstern  probabilities  in 
decision  theory  [-42}  and  the  concept  of  Turing  machines 
approximating  the  Universe  within  the  universal  artificial 
intelligence  model  of  Marcus  Hutter  [43} . 

We  now  briefly  discuss  some  of  the  quantum  phenom¬ 
ena  which  might  be  explained  by  the  CDFT  in  line  with 
a  deterministic  many-worlds  interpretation  of  quantum 
mechanics,  which  has  already  been  discussed  in  litera¬ 
ture  [44}.  Moreover,  we  propose  an  approach  to  demon¬ 
strate  unitary  evolution  for  quantum  fields. 

Heisenberg  uncertainty  principle.  There  is  a 
fundamental  limit  in  quantum  mechanics  to  determine 
both  position  and  momentum  of  a  particle  at  the  same 
time  [45[.  Werner  Heisenberg  first  introduced  this  prin¬ 
ciple  and  then  Earle  Kennard  proved  it  in  a  general  form 
of  <Jxap  >  -|,  where  h  is  the  Planck  constant  and  ax  and 
<jp  are  standard  deviations  of  position  and  momentum  of 
a  particle  [16}  . 

Howard  Robertson  proved  a  more  general  uncertainty 
relation  [47}  for  any  two  quantum  variables  A  and  B: 
&A&B  >  ±|<#A,BM|.  This  relation  shows  that  there 
is  a  fundamental  uncertainty  in  a  simultaneous  measure¬ 
ment  of  two  non-connnuting  quantum  variables.  This 
fact  is  as  a  manifestation  of  incomplete  description  of 
the  Universe  by  quantum  mechanics. 

The  many-worlds  interpretation  view  on  the  uncer¬ 
tainty  relations  is  that  every  quantum  measurement 
splits  the  Universe  so  that  observers  calculate  statistical 
correlations  of  two  consecutive  measurements  in  different 
splits  of  the  Universe.  The  simultaneous  measurement 
is  only  possible  for  commuting  quantum  operators  be¬ 
cause  they  have  common  eigenvectors  so  that  the  result 
of  the  measurement  for  both  operators  happens  in  the 
same  split  of  the  Multiverse. 

It  is  an  interesting  direction  for  future  research  for  the 
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CDFT  to  find  an  analog  of  uncertainty  relations  and  to 
check  them  computationally. 

Double-slit  experiment.  Classical  waves  can 
demonstrate  diffraction  and  interference,  but  quantum 
mechanics  has  a  unique  twist  on  top  of  that.  As  Richard 
Feynman  noted,  the  double-slit  experiment  demonstrates 
the  core  phenomenon  of  quantum  mechanics,  which  can¬ 
not  be  explained  in  any  classical  way  [48j.  In  this  ex¬ 
periment,  the  interference  is  destroyed  any  time  an  ex¬ 
perimenter  measures  which  path  is  taken  by  a  particle. 
David  Deutsch  argued  that  this  phenomenon  is  a  man¬ 
ifestation  of  the  many- worlds  interpretation  of  quantum 
mechanics  [H]. 

We  follow  the  same  explanation.  In  the  CDFT,  a 
quantum  measurement  happens  anytime  when  the  Mul- 
tiverse  splits  irreversibly.  We  conjecture  that  interfer¬ 
ence  happens  each  time  when  multiple  Universes  coincide 
everywhere  except  for  a  finite  causal  diamond  region  of 
D(v,w)  =  J+{v)  D  so  that  there  are  multiple 

possible  field  configurations  inside  D(v,w ),  so  that  the 
future  and  the  past  are  the  same. 

On  the  contrary,  the  interference  disappears  whenever 
splits  of  the  Multiverse  began  to  differ  forever  after  some 
point  along  the  world  line  of  an  observer.  The  role  of 
the  observer  is  passive.  She  notices  post-factum,  which 
particular  variant  of  the  Universe  has  happened. 

EPR  paradox  is  a  thought  experiment  that  first 
demonstrated  that  quantum  mechanics  is  not  a  complete 
description  of  physical  reality  In  this  experiment,  a 
measurement  of  a  spin  for  one  of  two  spatially  separated 
spin-singlet  particles  instantaneously  defines  the  spin  for 
the  other  particle.  Either  a  “spooky  action  at  a  distance” 
can  describe  this,  which  is  contrary  to  general  relativity, 
or  some  hidden  variables  define  full  deterministic  evolu¬ 
tion  of  each  particle  before  the  measurement. 

The  proposed  theory  has  to  preserve  perfect  correla¬ 
tions  between  the  fields  along  the  world  lines  of  particles 
of  the  EPR  pair  in  order  to  explain  the  paradox  in  a  local 
and  deterministic  way.  This  requirement  is  quite  strong, 
but  still  possible  in  principle.  There  must  be  a  split  of 
the  Multiverse  at  the  spacetime  event  when  the  two  par¬ 
ticles  separated  from  each  other.  The  two  variants  of  the 
Universe  must  have  identical  world  lines  of  two  particles 
except  for  their  spin  variables  until  the  next  split  of  the 
Multiverse  when  their  spin  projections  are  measured. 

Bell  theorem  provides  another  challenge  for  a  lo¬ 
cal  and  deterministic  theory.  The  original  Bell’s  work 
stated  that  any  local  hidden  variable  classical  model 
could  not  describe  all  phenomena  in  quantum  mechanics 
from  an  example  of  two  spin-^  particles  forming  an  EPR 
pair  (5lj] .  However,  later  it  was  noted  that  a  classical 
super-deterministic  model  could  allow  for  satisfaction  of 
Bell- type  experiments  [52j.  In  this  hypothetical  model, 
an  observer  itself  is  a  part  of  a  fully  deterministic  en¬ 
vironment,  so  her  decisions  to  choose  a  polarization  of  a 
measurement  device  depends  on  the  state  of  the  Universe. 
There  is  no  free  will  for  an  observer  in  such  interpretation 
of  the  paradox. 


A  three  particle  GHZ  experiment  provides  a  more 
clear  explanation  of  the  Bell  theorem  without  statisti¬ 
cal  inequality  relations  in  the  original  work.  It  can  be 
shown  that  a  quantum  state  of  three  spin-^  particles 
I  GHZ)  =  A=(|  tzUI  Ub\  t  z)c  -  I  U)a\  U)b\  iz)c) 
cannot  have  definite  values  for  x-  and  y-projections  of  all 
spins  simultaneously,  which  means  there  is  no  Universe 
where  hidden  variables  define  these  projections. 

In  Lev  Vaidman’s  interpretation  of  this  experi¬ 
ment  [5.3| ,  the  outcome  of  the  measurement  of  spins  of 
the  second  GHZ  quantum  state  along  the  x-direction  de¬ 
fines  the  choice  of  either  x-  or  y-direction  to  measure 
the  spin  of  the  first  GHZ  quantum  state.  Both  quantum 
mechanics  and  the  many-worlds  interpretation  allow  for 
only  16  possible  combinations  of  outcomes  for  six  spin- 
i  measurements  in  four  splits  of  the  Multiverse,  while  a 
free  choice  of  measurement  directions  for  the  first  GHZ 
quantum  state  leads  to  64  possible  combinations  of  the 
outcomes. 

This  interpretation  of  the  Bell  theorem  allows  keep¬ 
ing  a  local  hidden  variable  description  of  reality  without 
the  need  of  instant  and  non-local  “spooky  action  at  a  dis¬ 
tance”  .  However,  the  theory  still  has  to  have  non-locality 
in  terms  of  preserving  perfect  correlations  for  spatially 
separated  objects.  We  argue  that  the  CDFT  is  an  excel¬ 
lent  candidate  to  describe  physical  reality  in  a  local  and 
deterministic  way  leading  to  these  perfect  spatial  corre¬ 
lations  within  a  single  split  of  the  Multiverse. 

Quantum  entanglement  still  does  not  have  a 
good  explanation  beyond  the  Copenhagen  interpretation, 
which  allows  for  quantum  states  of  spatially  separated 
particles  that  can  be  described  by  a  joint  wavefunction, 
but  any  subset  of  these  particles  cannot  be  described  sep¬ 
arately  from  other  particles.  The  proposed  theory  pro¬ 
vides  an  interpretation  of  this  phenomenon  in  a  local  and 
deterministic  way. 

We  argue  that  quantum  entanglement  arises  from  non- 
computable  propagation  rules.  Without  them,  the  Mul¬ 
tiverse  would  have  all  possible  configurations  of  integer 
value  fields  and  directions  of  links  subject  to  a  unique  Big 
Bang  event  and  the  directed  acyclic  graph  property.  Non- 
computable  propagation  rules  set  additional  global  con¬ 
straints  that  limit  the  number  of  possible  local  splits  of 
the  Multiverse,  and  therefore  create  strong  long-distance 
correlations  between  discrete  fields  within  a  single  split 
of  the  Multiverse. 

Unitary  evolution  for  quantum  fields.  On  a  final 
note,  we  introduce  a  direction  for  future  research  toward 
a  quantum  held  theory  interpretation  of  the  CDFT  that 
will  have  unitary  evolution  for  quantum  fields.  We  use 
Gerard  ‘t  Hooft’s  ideas  for  a  deterministic  interpretation 
of  quantum  mechanics,  which  follows  from  the  cyclical 
evolution  of  discrete  systems  [Hi}. 

We  get  back  to  a  stationary  Universe  example  with  a 
nanotube  graph  to  demonstrate  unitary  evolution.  There 
are  microstates  |$j)  and  macrostates  |4/9)  for  each  charge 
value  q.  All  microstates  form  a  Markov  stochastic  matrix 
A4  for  probabilistic  transitions  from  one  microstate  to 


17 


another.  This  matrix  has  non-overlapping  probability 
paths  for  each  macrostate  with  charge  q. 

The  key  idea  is  to  find  a  unitary  matrix  U  in  the  space 
of  abstract  vectors  |$j),  that  has  the  same  transition 
probabilities  Pjk  =  Mkj  =  \Ukj\2- 

There  are  two  general  types  of  stochastic  matri¬ 
ces  that  have  this  property.  One  example  is  a 
Hadamard  matrix  in  N  dimensions:  Ujk  =  Hjk  = 
exp  [2ni(j  —  1  )(fc  —  1  )/N]/\/TSi.  This  Hadamard  matrix 
generates  probabilistic  transitions  from  each  state  equally 
likely  to  all  N  states.  Another  example  is  a  deterministic 
transition  matrix  that  has  all  zeros  and  a  single  one  for 
each  row  and  each  column.  The  proposed  theory  with 
the  postulated  principle  of  least  action  for  equiprobable 
sequential  propagation  has  elements  of  both  Hadamard 
type  evolution  and  deterministic  evolution. 

Here  is  an  example  of  how  a  particular  case  of  a  de¬ 
terministically  evolving  CDFT  model  can  have  unitary 
evolution  between  microstates  |$j).  We  consider  deter¬ 
ministic  sequential  growth,  so  that  there  is  an  algorithm 
that  takes  topology  of  a  microstate  1$,}  and  determin¬ 
istically  calculates  which  next  vertex  to  propagate.  If 
the  algorithm  is  complex  enough,  then  it  can  generate 
a  pseudo-random  sequence,  so  that  the  postulated  prin¬ 
ciple  of  least  action  still  holds  statistically.  This  deter¬ 
ministic  propagation  uniquely  defines  a  split  of  the  Mul- 
tiverse,  so  that  any  microstate  has  deterministic  future. 

Now  coming  back  to  the  stationary  Universe  exam¬ 
ple,  all  microstates  during  sequential  growth  will  have 
a  unique  immediate  descendant.  Some  microstates  will 
form  circles  of  consecutive  transitions,  and  some  mi¬ 
crostates  will  converge  to  these  circles  exactly  as  Gerard 
‘t  Hooft  proposed  in  his  work  on  the  deterministic  in¬ 
terpretation  of  quantum  mechanics  [54j.  From  the  anal¬ 
ogy  in  real  physics,  microstates  from  all  circles  form  true 
vacua,  and  the  remaining  microstates  form  false  vacua 
that  eventually  decay  to  one  of  the  true  vacua. 

Following  the  ideas  in  [Sdj.  circular  evolution  of  mi¬ 
crostates  in  a  certain  true  vacuum  is  described  by  an 
evolution  operator  with  a  unitary  matrix  Uop  that  has 
some  Hermitean  operator  with  a  Hermitean  matrix  Hop: 

Uop  =  e~iH°r5t  (32) 

That  is  how  a  Hilbert  space  for  the  quantum  mechan¬ 
ical  description  may  appear  in  the  proposed  theory.  Mi¬ 
crostates  |$i)  are  ontological  states  in  this  Hilbert  space 
that  allows  for  a  change  of  the  basis  to  any  linear  su¬ 
perposition  of  them.  This  approach  shows  the  way  how 
complex  numbers  appear  in  the  proposed  theory,  which 
is  essential  because  they  are  crucial  for  quantum  mechan¬ 
ical  description. 

A  static  Universe  with  deterministic  propagation  is  a 
minimal  case.  More  research  is  required  to  study  the 
relevance  of  general  CDFT  models  for  the  interpretation 
of  quantum  mechanics.  Two  key  questions  are  what  to  do 


with  an  expanding  Universe  and  whether  deterministic 
propagation  is  a  crucial  part  of  the  theory. 

V.  CONCLUSION 

We  propose  a  novel  approach  to  quantum  gravity  to 
unite  general  relativity  and  quantum  mechanics.  We 
start  with  a  core  postulate  that  a  complete  description 
of  physical  reality  has  to  be  causal,  discrete,  local,  and 
deterministic.  We  propose  to  consider  models  with  fun¬ 
damental  symmetries  of  isotropy,  CPT  invariance,  and 
charge  conservation.  The  resulting  theory  is  a  covariant 
combination  of  cellular  automata  and  the  causal  set  the¬ 
ory  on  self-similar  graphs.  The  beauty  of  this  theory  is 
that  it  is  static.  In  some  sense,  the  Universe  has  already 
happened,  and  time  is  just  an  illusion  that  follows  from 
causal  relations  for  the  fields  around  an  observer. 

We  studied  only  simple  computable  models  for  hun¬ 
dreds  of  Planck  epochs  from  the  Big  Bang  because  the 
computation  of  even  a  hundred  epochs  in  four  dimensions 
takes  weeks  on  a  single  computer.  The  theory  reproduces 
an  eternally  expanding  closed  Universe  with  a  physically 
relevant  scale  for  the  Hubble  constant  and  the  cosmolog¬ 
ical  term.  We  only  speculated  on  how  non-computable 
models  can  be  relevant  to  demonstrate  mass  and  inertia, 
and  to  reproduce  cosmological  parameters  for  an  early 
Universe  then  matter  played  a  dominant  role  in  space- 
time  expansion.  We  also  demonstrated  unitary  evolution 
only  for  a  stationary  Universe  with  deterministic  propa¬ 
gation.  These  are  two  major  week  areas  for  the  proposed 
theory  that  require  additional  research  in  the  future. 

Any  scientific  theory  has  to  be  falsifiable  and  to  have  an 
experimental  confirmation.  There  are  several  approaches 
to  that.  First,  one  may  try  to  explain  large-scale  phenom¬ 
ena,  e.g.,  the  patterns  of  discrete  fields  near  the  Big  Bang 
event  may  explain  the  evolution  of  cosmological  variables 
for  early  stages  of  expansion  of  the  Universe  [26j,  or  the 
underlying  structure  of  the  self-similar  graph  may  be  rel¬ 
evant  to  explain  anisotropy  and  angular  differences  in 
higher  moments  of  the  cosmic  microwave  background  ra¬ 
diation  [55] .  Second,  one  may  perform  an  ab  initio  calcu¬ 
lation  of  known  fundamental  constants,  e.g.,  by  defining 
energy  from  quantum  field  theory  approach  and  check¬ 
ing  its  correspondence  with  gravitational  energy  derived 
from  geometry.  Finally,  one  may  search  for  stable  pat¬ 
terns  that  may  represent  light  or  elementary  particles  and 
study  interactions  between  them  to  get  to  the  Standard 
Model  of  particle  physics  as  low  energy  and  long-distance 
approximation. 

There  is  still  much  room  to  prove  rigorously  that  the 
proposed  causal  discrete  field  theory  is  relevant  to  de¬ 
scribe  physical  reality.  If  this  is  a  right  approach  for 
quantum  gravity,  there  is  still  a  need  to  find  the  right 
combination  of  two  things:  an  underlying  self-similar 
graph  and  a  propagation  rule.  In  any  case,  the  proposed 
theory  opens  up  a  novel  way  to  describe  the  Universe  at 
all  scales. 
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